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(54) Road surface condition estimating apparatus and variation reduction processing apparatus 



(57) Disclosed is a road surface condition estimat- 
ing apparatus including; a converting portion (12, 14, 
16) having a table showing a relation among a vehicle 
speed (V). a braking force (Pc) and a resonance gain 
(Gd) for each road surface condition and converting a 
detected value of the input vehicle speed and braking 
force into a resonance gain (Gd1, Gd2, Gd3) on the 
basis of the table, a subtracting device (18, 20, 22) for 
computing each of differences between the resonance 
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gain (Gd-Gd1, ...) for each road surface condition con- 
verted by each of the converting portions and the 
detected value of the input resonance gain (Gd). and a 
minimum value selecting portion (24) selecting the min- 
imum value among the computed values in the subtract- 
ing device (18, 20, 22) and outputting an information of 
the road surface condition corresponding to the convert- 
ing portion (12, 14, 16) giving the minimum value. 
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Description 

BACKGROUND OF THE INVENTION 
5 Held of the Invention 

The present invention relates to a road surface condition estimating apparatus lor estimating the condition of a road 
surface on which a vehicle is running, and more specifically to a road surface condition estimating apparatus which is 
provided with a correlation among a plurality of wheel motion amounts expressing a road surface \i characteristic for 
io each road surface condition, and which estimates and outputs the road surface condition on the basis of the correlation. 
The present invention also relates to a variation reduction processing apparatus and more specifically to a variation 
reduction processing apparatus which can be applied to a road surface condition estimating apparatus, an ABS (anti- 
lock braking system) apparatus, a driving force controlling apparatus, and the like, and to a road surface condition esti- 
mating apparatus to which the variation reduction processing apparatus is applied. 

75 

Description of the Related Art 

Conventionally, an anti-lock brake control device has been proposed wherein, immediately prior to a vehicle's 
wheel shifting to a locked up state when the coefficient of friction \i between the wheel and the road surface surpasses 
20 a peak value of ji, (hereinafter referred to as "a peak n" or "a peak \x value") this device can prevent the vehicle's wheels 
from being locked up by reducing the braking torque acting on the wheel and by controlling the coefficient of friction u 
to follow this peak value. 

Meanwhile, when a vehicle is running at a certain speed and the brakes are applied, then the vehicle's wheels slip 
on the road surface, ft is known that the coefficient of friction ji between the wheels and the road surface varies, as is 
25 shown in Fig. 21 in relation to the slip ratio X expressed in formula (1) below. 

* = (V-V W )/V (1) 

in which V is vehicle speed (conversion into an angular speed) and V w is a wheel speed, and accordingly (V - 
30 V w ) corresponds to a slip speed Av. 

As shown in Fig. 21 , in this \i-X characteristic, the coefficient of friction ji reaches the peak value at a certain slip 
speed (the area A2 in Fig. 21). 

In Japanese Patent Application Laid-open (JP-A) No. 1-249559, there has been proposed the following anti-lock 
brake controlling apparatus. The anti-lock brake controlling apparatus computes the slip ratio from the approximate 
35 vehicle speed and the detected wheel speed using the formula (1). The anti-lock brake controlling apparatus controls 
the braking force in such a manner that the computed slip ratio substantially coincides with a reference slip ratio (a slip 
ratio giving the peak y) which has been previously set at a fixed value. Accordingly, the coefficient of friction \i is contin- 
ually controlled at the peak \i value. 

However, the slip ratio at which the coefficient of friction \i reaches the peak u value varies in accordance with the 
40 road surface condition (the type of road surface on which the vehicle is running and the condition thereof). Accordingly, 
when performing follow-up control according to the fixed reference slip ratio as in the prior art described in the publica- 
tion mentioned above, the stopping distance may become too long depending on the road surface, or else the brakes 
may be controlled so that the coefficient of friction \x goes over the peak \i value, whereby there is the concern that the 
wheels may lock. 

45 In order to solve this drawback, it is necessary to estimate the road surface condition and change the reference slip 
ratio in accordance with the estimated road surface condition. However, in the conventional ait, there is no technique 
for accurately estimating the road surface condition. 

SUMMARY OF THE INVENTION 

50 

In view of the facts mentioned above, an object of the present invention is to provide a road surface condition esti- 
mating apparatus for accurately estimating a road surface condition. 

In order to achieve the object mentioned above, in accordance with a first invention of the present invention, there 
is provided a road surface condition estimating apparatus provided with a correlation among a plurality of wheel motion 
55 amounts expressing the characteristic of a coefficient of friction n between a wheel and a road surface with respect to 
a slip speed for each road surface condition, comprising: 

memory means for storing a gradient of the coefficient of friction \x with respect to the slip speed or a physical 
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amount relating to said gradient as one of the wheel motion amounts expressing said correlation; and 
road surface condition estimating means for estimating and outputting a road surface condition on the basis of com- 
parison between a converted value of the wheel motion amount for each road surface condition obtained by 
respectively converting at least one of detected values of the plurality of wheel motion amounts in accordance with 
5 said stored correlation and a detected value of the same amount as said wheel motion amount. 

In this case, a characteristic of the coefficient of friction \i between the wheel and the road surface with respect to 
the slip speed Av corresponds to a Av - n characteristic expressed in Fig. 17, for example, and includes all the other 
equivalent characteristics (hereinafter, referred to as "road surface ji characteristic"). As the characteristic equivalent to 

w the Av - fi characteristic, for example, there are a characteristic shown in Fig. 21 and expressed by a slip ratio in place 
of the slip speed, a characteristic expressed by a braking force relating to \i foiW: W is a wheel bad) and a braking 
torque QiWr: r is an effective radius of the wheel) in place of the coefficient of friction \i. and the like. The road surface 
li characteristic varies corresponding to the road surface condition as shown in Fig. 2A. 

Further, the wheel motion amount corresponds to a physical amount relating to the wheel motion, which is neces- 

15 sary for expressing the road surface u characteristic, for example, there are vehicle speed, a wheel speed, a wheel 
deceleration (or acceleration), a slip speed, a braking force, a braking torque, a braking pressure (a wheel cylinder pres- 
sure) and the like as well as the physical amounts expressing slipperiness of the wheel such as a gradient of the road 
surface ji with respect to the slip ratio, a gradient of the braking torque with respect to the slip speed, and the like. 
These wheel motion amounts are selected in such a manner as to necessarily and sufficiently express the road sur- 

20 face ji characteristic respectively becoming a fixed relation for each road surface condition, and accordingly, the corre- 
lation between the wheel motion amounts directly or indirectly expresses the road surface \i characteristic for each road 
surface condition. 

For example, the road surface \i characteristic in Fig. 2A can be at first expressed for each road surface condition 
by the correlation between the braking force P c as another wheel motion amount and the slip ratio X, as shown in Fig. 
25 2B. Then, the road surface ii characteristic in Fig. 2A can be necessarily and sufficiently expressed by the braking force 
P c , the vehicle speed and the gradient of the road surface ii as a three-dimensional characteristic for each road sur- 
face condition, as shown in Fig. 2C, considering the fact that the gradient Gd of the road surface ii with respect to the 
slip ratio X, as well as the \i value, has a peculiar value for each road surface condition in Fig. 2A, and using the relation 
in Fig. 2B. 

30 In this case, since the slip ratio X can be expressed by the vehicle speed and the wheel speed, in the case of using 
the braking force as one of the wheel motion amounts, another wheel motion amount is required, so that in an embod- 
iment shown in Fig. 2C, the vehicle speed is employed. 

In Fig. 2C, the three-dimensional characteristic is expressed by the relation between the braking force and the road 
surface n gradient in the case of each of the vehicle speed V of 20 km/h, 40 km/h and 60 km/h. From this drawing, 

35 it is understood that the relation among these three wheel motion amounts is different in accordance with the road sur- 
face condition. 

In accordance with the present invention, the correlation between a plurality of wheel motion amounts expressing 
these road surface ii characteristics is provided for each road surface condition. In this case, the correlation can be 
expressed, for example, by a conversion table between the input and output data or the like, and in the case of the char- 
40 acteristic of Fig. 2C, the table is previously formed on the basis of the data obtained by detecting the relation between 
the braking force and the road surface ii gradient for each vehicle speed and each road surface condition. 

Here, in the case of estimating the road surface condition when the wheel motion amount such as the vehicle speed 
and the like is constant, it is possible to reduce a number of the wheel motion amounts expressing the correlation. 
In this case, the present invention comprises detecting means for detecting a plurality of wheel motion amounts 
45 expressing a characteristic of a coefficient of friction \i between a wheel and a road surface with respect to a slip speed ; 

converting means for converting a detected value of at least one of the plurality of wheel motion amounts into a con- 
verted value of another of the plurality of wheel motion amounts for each of a plurality of road surface conditions on 
the basis of a correlation between the plurality of wheel motion amounts for each of a plurality of road surface con- 
so ditions provided in said converting means; and 

road surface condition estimating means for estimating a road surface condition by comparing said converted value 
with a detected value of said another of the plurality of wheel motion amounts. 

Accordingly, the road surface condition estimating means in accordance with the present invention can be struc- 
55 tured in such a manner that a plurality of wheel motion amounts can be input. For example, as shown in Fig. 1 , the struc- 
ture is made such that the vehicle speed V detected by the vehicle speed detecting means, the gradient of the road 
surface ii detected by the road surface \i gradient detecting means and the braking force P c detected by the braking 
force detecting means can be input to the road surface condition estimating means. 
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When the detected values of the plurality of wheel motion amounts are input, the converting means converts at 
least one (one or a plural number) of the detected values of the plurality of input wheel motion amounts on the basis of 
the correlation for each of the road surface conditions. As a result, the converted values of the wheel motion amounts 
other than the wheel motion amount which is a subject of conversion can be obtained for each road surface condition. 

5 For example, in the structure shown in Fig. 1 , in the case of converting the vehicle speed and the braking force in 
accordance with the correlation shown in Fig. 2C, the converted values of the gradient of the road surface \i correspond- 
ing to the detected value of the vehicle speed V = 20 km/h and the detected value Pco of the braking force respectively 
become G^ , G^ and for a dry road surface, a snow covered road surface and a ice covered road surface. 

Then, the road surface condition estimating means estimates the road surface condition on the basis of the com- 

io parison between the converted values of the other wheel motion amounts for each of the road surface condition 
obtained in this manner and the corresponding detected values of the same wheel motion amounts. For example, in the 
case of the embodiment mentioned above, the converted values Qdi. Gd2 and 0^ of the gradient of the road surface 
u for each of the road surface conditions are compared with the detected value of the gradient of the road surface 
|i detected by the road surface \i gradient detecting means and the road surface condition is estimated and output from 

is the result of the corrparison. 

In this comparison, for example, a difference between the detected value Gdo and each of the converted values 
Gdi . Gd2 and Gd3 is computed, and the road surface condition giving the smallest difference (absolute value) is output. 
For example, in the case the difference between G^ and is the smallest, it estimates that the road surface condition 
is "snow covered road surface". 

20 Further, the apparatus is characterized in that one of a plurality of wheel motion amounts expressing the correlation 
is a gradient of the coefficient of friction \i with respect to the slip speed or a physical amount relating to the gradient. 

Since the gradient of the coefficient of friction \i with respect to the slip speed (hereinafter referred to as The gradi- 
ent of the road surface \i w or The road surface \x gracfienT) has a peculiar value in correspondence to the road surface 
condition as already mentioned above, it becomes an effective wheel motion amount at a time of estimating the road 

25 surface condition. Further, it is hard to directly detect the coefficient of friction u, while an effective method of computing 
by using an online system identification method from time series data of the wheel speed (the invention described in 
Japanese Patent Application No. 8-21 8828) can be applied to obtain the gradient of the road surface u. Accordingly, the 
road surface condition can be more accurately computed by using the gradient of the road surface \i. 

In this case, as a physical amount relating to a gradient of the coefficient of friction \i with respect to the slip speed, 

30 for example, there are a gradient of a braking torque (or a braking force) with respect to the slip speed (or the slip ratio) 
and the like. 

As mentioned above, in accordance with the present invention, an accurate road surface condition can be esti- 
mated on the basis of the correlation between a plurality of wheel motion amounts expressing a characteristic of the 
coefficient of friction n between the wheel and the road surface with respect to the slip speed and the detected values 

35 of a plurality of wheel motion amounts, and a vehicle stabilizing control can be more accurately performed on the basis 
of the estimated road surface condition. 

Here, in accordance with the present invention, the structure may further comprise micro vibration means for 
minutely vibrating a braking force by a resonance frequency of a vibration system constituted by a vehicle body, a wheel 
and a road surface, in which a physical amount relating to the gradient of the coefficient of friction \i may be set to a 

40 resonance gain, that is, a ratio of the micro amplitude of a resonance frequency component of the wheel speed with 
respect to the micro amplitude of the braking force at a time of minutely vibrating the braking force by means of the 
micro vibration means. 

Accordingly, a vibration phenomenon in the wheel at a time when the vehicle comprising a vehicle body of a weight 
W v is running at vehicle speed that is, a vibration phenomenon in a vibration system constituted by the vehicle body, 
45 the wheel and the road surface is considered by referring a model shown in Fig. 16 which is equivalently modeled by a 
wheel rotation axis. 

In the model shown in Rg. 16, the braking force acts on the road surface through a surface of a tread 115 of a tire 
being in contact with the road surface. Here, since the braking force actually acts on the vehicle body as a reaction (a 
braking force) from the road surface, an equivalent model 1 17 of the vehicle body weight W v in accordance with the 
so rotation axis conversion becomes connected to an opposite side of the wheel 1 13 through a friction element 116 (the 
road surface u) between the tread of the tire and the road surface. This is similar to the case that the weight of the vehi- 
cle body can be modeled by a great inertia under the wheel, that is, a mass on the opposite side of the wheel, as in a 
chassis dynamo apparatus. 

In Rg. 1 6, when an inertia of the wheel 1 1 3 including a tire rim is set to be J w a spring constant of a spring element 
55 114 between the rim and the tread 1 15 is set to be K, a radius of the wheel is set to be R, an inertia of the tread 115 is 
set to be J t , a coefficient of friction of a friction element 116 between the tread 1 15 and the road surface is set to be \i. 
and an inertia of the equivalent model 1 1 7 of the vehicle body weight in accordance with the rotation axis conversion is 
set to be a transmission characteristic from a braking torque T b ' generated by the wheel cylinder pressure to the 
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wheel speed is computed by a formula of the wheel motion and is expressed as follows. 

^oJJ<fW^J^y+a(J t +J Y )R 2 Wy+JJCs+aKR 2 W v } T b 
<O w = 5 5 5 (2) 



In which reference symbol s is an operator for Laplace transform. 

When the tire grips on the road surface, in the case that the tread 1 15 and the vehicle equivalent model 117 are 
10 considered to be cfirectly connected to each other, an inertia of the sum of the vehicle equivalent model 117 and the 
tread 1 15 and an inertia of the wheel 113 resonate. Accordingly, the vibration system can be considered to be a wheel 
resonance system constituted by the wheel, the vehicle body and the road surface. The resonance frequency ©oo of the 
wheel resonance system at this time is expressed in the transmission characteristic of the formula (2) as follows. 



15 



<*» = V«J w -M t v)K/J W (J v)}** H> 



The state corresponds to an area A1 prior to reaching the peak \i in Fig. 17 (Fig. 21). 

On the contrary, in the case that the coefficient of friction ji of the tire approaches the peak u, the coefficient of fric- 
20 tion ^ of the tire surface is hard to be changed with respect to the slip ratio, and the component in accordance with the 
vibration of the inertia of the tread 1 15 does not affect the vehicle body equivalent model 117. Accordingly, the tread 
1 1 5 and the vehicle equivalent model 1 1 7 are equivalents separated, and the tread 1 1 5 and the wheel 1 1 3 generate 
resonance. The wheel resonance system at this time can be considered to be constituted by the wheel and the road 
surface, and the resonance frequency ©oo at this time equals to a value obtained by setting the vehicle body equivalent 
25 inertia Jy to be 0 in the formula (3). Accordingly, the value &x> is expressed by the following formula. 



= VttJw-KJtJK/J^,)}^ (4) 

30 The state corresponds to an area A2 near the peak \i in Fig. 17 (Fig. 21). 

Comparing the formula (3) with the formula (4), when it is supposed that the vehicle body equivalent inertia Jy is 
larger than the wheel inertia J w and the tread inertia Jt, the resonance frequency ©oo' of the wheel resonance system in 
the case of the formula (4) is shifted to the higher frequency side than <d°o in the formula (3). Accordingly, it is possible 
to determine the friction state between the wheel and the road surface on the basis of the physical amount reflecting 
35 the change of the resonance frequency in the wheel resonance system. 

Then, in accordance with the present invention, as the physical amount for reflecting the change of the resonance 
frequency of this kind, a resonance gain Qd as that mentioned below is introduced as a wheel motion amount. 

At first, when the micro vforation means minutely vibrates the braking force at the resonance frequency cdoo (the for- 
mula (3)) in the vibration system comprising a wheel, a vehicle body and a road surface (in this case, the braking pres- 
40 sure P b is supposed to be minutely vibrated), the wheel speed ©* is also minutely vibrated at the resonance frequency 
©ao around the average wheel speed. Here, in the case of setting the micro amplitude of the braking pressure P b at the 
resonance frequency ©°o at this time to be P v and setting the micro amplitude of the wheel speed at the resonance fre- 
quency ©°o to be a>w the resonance gain Qd is expressed by the following formula. 

45 G^ffl^y (5) 

In this case, the resonance gain Qd can be considered to be a vibration component of the resonance frequency ©oo 
in a ratio {(oJP b ) of the wheel speed <a* with respect to the braking pressure and expressed by the following for- 
50 mula. 

G d = ((<D v /P b )|s = j(Doo) (6) 

55 Since the resonance gain Qd corresponds to a vibration component of the resonance frequency a>oo in the ratio 
(©w/Pb) as shown in the formula (6), the resonance frequency is shifted to the value of ©oo' when the friction state 
reaches the area adjacent to the peak u, so that it is suddenly reduced. Accordingly, the resonance gain Qd corre- 
sponds to the wheel motion amount sensitively reflecting the friction state defining the characteristic of the road surface 
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Next, explanation will be given to the fact that the resonance gain is a physical amount equivalent to a braking 
torque gradient (a gradient of the braking torque with respect to the slip speed) which corresponds to one expression 
of the gradients of the road surface \i. 
s As is already mentioned, the characteristic in Fig. 1 7 shows that a functional relation in which the coefficient of fric- 
tion *i reaches the peak at a certain slip speed can be established between the slip speed Av and the coefficient of fric- 
tion ii between the wheel and the road surface. 

In this case, since the wheel speed is minutely vibrated as a result of minutely vibrating the braking pressure by 
means of the micro vibration means, the slip speed also minutely vibrates around a certain slip speed. Here, in the road 
10 surface condition having the characteristic of Fig. 17, a change of the coefficient of friction \x with respect to the slip 
speed Av at a time of minutely vibrating the slip speed around a certain slip speed will be considered. 

At this time, the coefficient of friction n on the road surface can be approximately expressed by the following formula 

W- 

15 \i = n 0 +aRAv (7) 



Since the change of the slip speed due to the micro vibration is small, it can be approximately expressed by a 
straight line having a gradient aR. 
20 In this case, with substituting the formula (7) for a braking torque T b = fiWR generated by the coefficient of friction 
ji between the tire and the road surface, the following formula (8) can be obtained. 



25 



30 



T b = jiWR = ji 0 WR+aR 2 AvW (8) 

In which reference symbol W is a wheel load. When both sides in the formula (8) are differentiated by Av for one 
stage, the following formula (9) can be obtained. 

dT h 2 

Wv =aRW < 9 > 



Accordingly, it is shown by the formula (9) that the braking torque gracBent (dTt/Av) equals to aR^W. 

On the contrary, since the braking torque T b ' is in proportion to the braking pressure Pb, the resonance gain Gd is 
35 in proportion to a vibration component of the resonance frequency <doo in the ratio of the wheel speed cow with respect 
to the braking torque T b ' (o\ v /T b '). Accordingly, the resonance gain can be expressed by the following formulas (10), 
(1 1) and (1 2) in accordance with the transmission characteristic of the formula (2). 



40 



45 



» jA + ctB (j:an imaginary unit) 



(10) 



50 

in which, 

J A - Jf-tJy+J Wt J B = t/ f +Jy (11) 

J A J V J W K J A J V J W K 
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Generally, in the formula (12), since the formula (13) can be established, the formula (14) can be obtained in 
accordance with the formulas (9) and (10). 

| A | = 0.012 « | B | = 0.1 (13) 

5 

Wv"°* < 14 > 

10 Accordingly, the gradient of the braking torque T b with respect to the slip speed Av is in proportion to the resonance 
gain Qd. 

As mentioned above, it is shown that the resonance gain Qd is a physical amount equivalent to the gradient of the 
road surface \i (the braking torque gradient), and it is understood that the road surface condition can be accurately esti- 
mated on the basis of the resonance gain G^. Further, since the resonance gain Qd is a wheel motion amount which 
15 accurately expresses the friction state on the basis of the resonance characteristic of the vibration system, the present 
invention can more accurately estimate the road surface condition. 

(The Other Aspects of the Present Invention) 

20 In accordance with a first aspect of the first invention of the present invention, the structure may be made such that 
said detecting means detects a wheel speed as one of said plurality of wheel motion amounts, 

said con-elation includes a correlation defined for each vehicle speed and a correlation between two wheel motion 
amounts other than the vehicle speed, 
25 said converting means converts a detected value of any one of said two wheel motion amounts into a converted 
value of the other wheel motion amount for each of a plurality of slip speeds on the basis of said correlation 
between said two wheel motion amounts, 

said slip speed calculating means calculates a slip speed on the basis of said converted value and a detected value 
of said the other wheel motion amount, 
30 said vehicle speed calculating means calculates the vehicle speed on the basis of the slip speed calculated by said 
slip speed calculating means and the wheel speed detected by said detecting means, and 
said road surface condition estimating means estimates the road surface condition on the basis of the vehicle 
speed calculated by said vehicle speed calculating means and said correlation defined for each vehicle speed. 

35 In this case, a second table showing a correlation between two wheel motion amounts other than the vehicle speed 
is provided for each slip speed, the slip speed calculating means calculates the slip speed on the basis of comparison 
between the converted value of the other wheel motion amount for each slip speed obtained by converting one detected 
value of the two wheel motion amounts for each slip speed in accordance with the second table and the detected value 
of the other wheel motion amount, and the vehicle speed calculating means calculates the vehicle speed on the basis 

40 of the slip speed calculated by the slip speed calculating means and the detected value of the input wheel speed. In this 
case, one of the plurality of wheel motion amounts showing the correlation in accordance with the table corresponds to 
the vehicle speed, and the road surface condition estimating means estimates the road surface condtion by using the 
vehicle speed calculated by the vehicle speed calculating means. 

The first aspect is provided with the second table showing the correlation between two wheel motion amounts other 

45 than the vehicle speed for each slip speed. The second table can be made by a procedure mentioned below in the case 
that the two wheel motion amounts other than the vehicle speed are set to be the braking force P c and the gradient Qj 
of the road surface \i. 

At first, in each of the cases that the road surface condition and the vehicle speed (the vehicle speed) as parame- 
ters are variously changed, the change of the slip speed Av with respect to the braking force P c is computed, and then 
so the relation shown in Rg. 4 can be obtained. Next, in each of the cases that the road surface condition and the vehicle 
speed as parameters are variously changed, the change of the gradient of the road surface \i with respect to the slip 
speed Av is calculated, and then the relations shown in Fig. 5A (Dry) and Rg. 5B (Snow) can be obtained. 

Further, in accordance with the relations in Figs. 4, 5A and 5B, P c and Gd with respect to the same slip speed Av 
are calculated for each of the vehicle speeds, and then the relations shown in Rgs. 6A and 6B can be obtained for each 
55 of the road surface conditions. The relation between P c and Qd for each of the slip speeds shown by Rgs. 6A and 6B 
forms the second table. In this case, a line of each of the slip speeds Av in Rgs. 6A and 6B is fixed without relation to 
the vehicle speed and the road surface. 

Further, on the basis of Rgs. 6A and 6B, for each of the vehicle speeds (20 knVh, 40 km/h and 60 km/h), the relation 
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between P c and as shown in Figs. 7A to 7C is obtained for each road surface condition. The relation in Figs. 7A to 
7C corresponds to the table for each road surface condition in the case of setting a plurality of wheel motion amounts 
to be the vehicle speeds V, P c and G#. 

In the first aspect, for example, in accordance with the second table, the slip speed is calculated on the basis of the 

5 comparison between the converted value of the gradient of the road surface u for each slip speed obtained by con- 
verting the detected value of the braking force P c for each slip speed and the detected value of the gradient Gd of the 
road surface p. In other words, it is calculated which line of the slip speed Av in Figs. 6A or 6B is closest to the relation 
between the detected P c and Gd, and the Av value of the closest line is taken as the slip speed. 

Since the slip speed can be calculated in the above manner, the vehicle speed calculating means in accordance 

10 with the first aspect calculates the vehicle speed on the basis of the computed slip speed and the detected value of the 
input wheel speed. And the road surface condition estimating means estimates the road surface corxfition on the basis 
of the table by using the vehicle speed computed by the vehicle speed calculating means. In the case of the table shown 
in Figs. 7A to 7C, it is calculated which relation of the road surface condition is closest to the relation between the 
detected P c and Gd among the relations between P c and Qd in the table shown in Rgs. 7A to 7C corresponding to the 

is computed vehicle speed, and the closest relation is taken as the road surface condition. 

As mentioned above, in accordance with the first aspect, in spite of having the table including the vehicle speed as 
one of the wheel motion amounts, it is not necessary to detect the vehicle speed, and in place thereof, the wheel speed 
is detected. Accordingly, the road surface condition estimating means in accordance with the first aspect can be struc- 
tured, for exanrple, as shown in Fig. 3, such that the wheel speed V w detected by the wheel speed detecting means, 

20 the gradient of the road surface \i detected by the road surface \i gradient detecting means and the braking force P c 
detected by the braking force detecting means can be input to the road surface condition estimating means. 

Since it is hard to accurately estimate the vehicle speed in comparison with the wheel speed, it is possible to more 
accurately estimate the road surface condition by estimating the vehicle speed on the basis of the second table and the 
wheel speed as in the first aspect, in comparison with the case of Fig. 1 in which the vehicle speed is input as the 

25 detected value. 

Further, in accordance with the first aspect, since the accurate slip speed can be obtained, it is possible to more 
accurately control by applying the invention to an anti-lock brake control apparatus for controlling the sip ratio to coin- 
cide with the reference slip ratio. 

In each of the inventions mentioned above, the road surface condition is estimated. In addition, the structure may 
30 be made to further comprise peak value searching means for searching a peak \i value corresponding to the road sur- 
face condition estimated by the road surface condition estimating means from a peak u value (u max value) for each 
road surface previously stored as data. 

Accordingly, since the road surface condition and the peak u value can be recognized, the limit value of the braking 
force (driving force) can be clarified, so that it is possible to change the control gain and to set a control target value in 
35 a vehicle stabilizing control such as VSC (vehicle stability control), ABS, TRC (traction control) and the like, and to give 
an alarm for the road surface condition to the driver. 

Further, the structure may be made such that a maximum braking force is provided for each of a plurality of road 
surface conditions and a plurality of vehicle speeds, and 

40 said detecting means detects at least vehicle speed and a braking force as said plurality of wheel motion amounts, 
and wherein said road surface condition estimating apparatus further comprises: 

maximum braking force selecting means for selecting the maximum braking force corresponding to said estimated 
road surface condition and said detected vehicle speed on the basis of the maximum braking force provided for 
each of said plurality of road surface conditions and said plurality of vehicle speeds, the road surface condition esti- 
45 mated by said road surface condition estimating means and the vehicle speed detected by said detecting means, 
and 

subtracting means for calculating the difference between the maximum braking force selected by said maximum 
braking force selecting means and the braking force detected by said detecting means. 

so In other words, the road surface condition estimating apparatus may further comprise braking force surplus calcu- 
lating means for computing the maximum braking force corresponding to the vehicle speed and the road surface con- 
dition estimated by the road surface condition estimating means from the maximum braking force for each of the road 
surface conditions and the vehicle speeds previously stored as data and for calculating and outputting the difference 
between the maximum braking force and the input value of the detected braking force (a second aspect of the present 

55 invention). 

Here, the invention in accordance with the second aspect can be structured, for example, as shown in Fig. 8. As 
shown in the drawing, the vehicle speed V, the gradient of the road surface u and the braking force P c are input to 
the road surface condition estimating means, and as already mentioned, the road surface condition estimating means 
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estimates the road surface condition on the basis of the wheel motion amount and the table. Then, it transmits the road 
surface condition to the braking force surplus calculating means of the present invention. 

The braking force surplus calculating means comprises maximum braking force selecting means for selecting the 
maximum braking force corresponding to the road surface condition estimated by the road surface condition estimating 
5 means and the vehicle speed detected by the vehicle speed detecting means from the maximum braking force for each 
of the road surface conditions and the vehicle speeds previously stored as the data, and subtracting means for subtract- 
ing the braking force detected by the braking force detecting means from the calculated maximum braking force. 

In this case, the maximum braking force has an inherent value for each road surface condition, and the inherent 
value is different for each vehicle speed. For example, the maximum braking force (P c Max1, PcMax2 and P c Max3) 
w changes in correspondence to the road surface condition (Dry, Snow and Ice) and the vehicle speed (40 km/h and 60 
km/h), as shown in Figs. 9A and 9B. The maximum braking force selecting means previously stores the maximum brak- 
ing force for each of the road surface conditions and the vehicle speeds as shown in Figs. 9A and 9B as the data, and 
therefore can select the maximum braking force correspond ng to the input road surface condition and the vehicle 
speed. Then, the detected braking force is subtracted from the calculated maximum braking force by the subtracting 
is means, so that a surplus degree of the braking force can be computed. The surplus degree is an index of a degree of 
the braking force which can be applied to the wheel. 

Here, in the case of applying the second aspect to the invention of the first aspect, it is possible to more accurately 
calculate the maximum braking force by using the vehicle speed calculated by the vehicle speed calculating means. 

As mentioned above, in accordance with the second aspect of the present invention, since the surplus degree of 
20 the braking force can be known, it ts possible to stably run the vehicle by controlling in such a manner as to preferentially 
divide the braking force and the driving torque to the wheel having a high surplus degree. 

In order to achieve the object mentioned above, a second invention comprises first detecting means for detecting 
a first physical amount expressing a braking condition of a wheel of a vehicle; 

25 second detecting means for detecting a vehicle speed of said vehicle; 

third detecting means for detecting a second physical amount which represents slipperiness of said wheel and 
which is determined in accordance with a road surface condition on which said wheel is running, said first physical 
amount, and said vehicle speed; 

calculating means for calculating said second physical amount for each of a plurality of predetermined road surface 
30 conditions on the basts of said first physical amount and said vehicle speed; 

comparing means for comparing a detected track of said second physical amount detected by said third detecting 
means with respect to said first physical amount with each of calculated tracks of said second physical amounts 
calculated for each of the plurality of road surface conditions by said calculating means with respect to said first 
physical amount; and 

35 estimating means for estimating a road surface condition corresponding to a track closest to said detected track 
among said calculated tracks as an actual road surface condition in which said wheel is running on the basis of the 
comparison results of said comparing means. 

The first detecting means detects the first physical amount expressing the braking condition for the wheel of the 
40 vehicle. Here, the first physical amounts include the braking force to the wheel and the slip speed. 
The second detecting means detects the speed of the vehicle (the vehicle speed). 

The third detecting means detects the second physical amount expressing the slipperiness of the wheel. The sec- 
ond physical amount is defined in correspondence to the road surface condition in which the wheel is running, the first 
physical amount and the speed of the vehicle. Concretely speaking, the second physical amount is a gradient of a coef- 

45 ficient of friction ji between the road surface and the wheel or of a braking force or the braking torque against the wheel 
with respect to the slip speed. Accordingly, when the gradients are steep, the peripheral speed of the wheel is high and 
the wheel is hard to slip. On the contrary, when the gradients are gentle, the peripheral speed of the wheel is low, and 
the wheel is easy to slip. Accordingly, these gradients express the slipperiness of the wheel. 

The calculating means calculates the second physical amount for each of a plurality of predetermined road surface 

so conditions on the basis of the first physical amount detected by the first detecting means and the vehicle speed detected 
by the second detecting means. The calculating means may previously calculate the relation defined in correspondence 
to the vehicle speed between the first physical amount and the third physical amount for computing the second physical 
amount tor each of a plurality of road surface conditions, may calculate the third physical amount for each of a plurality 
of road surface conditions on the basis of the detected first physical amount and the detected vehicle speed, and may 

55 calculate the second physical amount for each of a plurality of road surface conditions from the calculated third physical 
amount Further, the calculating means may previously calculate the relation defined in correspondence to the vehicle 
speed between the first physical amount and the second physical amount for each of a plurality of road surface condi- 
tions, and may calculate the second physical amount for each of a plurality of road surface conditions on the basis of 
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the detected first physical amount and the vehicle speed. 

The comparing means compares the track of the second physical amount detected by the third detecting means 
with respect to the first physical amount expressing the braking condition of the wheel with each of tracks of the second 
physical amounts calculated for each of a plurality of road surface conditions by the calculating means with respect to 

5 the first physical amount. The estimating means estimates the road surface condition corresponding to the track closest 
to the detected track among the calculated tracks as the actual road surface condition in which the wheel is running on 
the basis of the comparison result of the comparing means. 

As mentioned above, the second physical amount is detected and the track of the detected second physical amount 
is compared with each of the tracks of the second physical amount calculated in correspondence to each of a plurality 

10 of road surface conditions. Accordingly, the road surface condition corresponding to the selected track corresponds to 
the actual road surface condition in which the wheel is running when selecting the track closest to the detected track 
among the tracks calculated in correspondence to each of a plurality of road surface conditions on the basis of the com- 
parison result. Thus, that it is possible to accurately estimate the actual road surface condition. 

In this case, the second physical amount mentioned above includes an error due to various causes. 

15 Accordingly, said first detecting means, said second detecting means and said third detecting means respectively 
detect the first physical amount, the vehicle speed, and the second physical amount at each predetermined time from 
the time braking is started on a wheel, 

said calculating means calculates the second physical amount for each of a plurality of road surface conditions at 
20 each time the first physical amount changes by a predetermined amount from the time braking is started on the 
wheel, and 

said comparing means compares the second physical amount detected at said predetermined time, with each of 
the second physical amounts calculated for each of the plurality of road surface conditions and at each time the first 
physical amount changes by a predetermined amount, at each time the first physical amount changes by a prede- 
25 ter mined amount from the time braking is started on a wheel. 

As a result, said comparing means compares the second physical amount detected at said predetermined time, 
with each of the second physical amounts calculated for each of the plurality of road surface conditions and at each time 
the first physical amount changes by a predetermined amount, at each time the first physical amount changes by a pre- 
30 determined amount from the time braking is started on a wheel, it is possible to reduce an influence due to the error, so 
that it is possible to more accurately estimate the road surface condition. 

In order to achieve the object mentioned above, a third invention comprises first detecting means for detecting a 
first physical amount expressing a braking condition or a driving condition of a wheel of a vehicle; 

35 second detecting means for detecting a second physical amount which represents sfipperiness of said wheel and 
varies in accordance with a change in said first physical amount; and 

reduction processing means for reducing variations in the second physical amount in accordance with said change 
in the first physical amount. 

40 The first detecting means detects the first physical amount expressing the braking condition or the driving condition 
of the wheel in the vehicle. In this case, the first physical amount includes the braking force against the wheel and the 
slip speed. 

The second detecting means detects the second physical amount expressing the slipperiness of the wheel. The 
second physical amount is concretely a gradient of a coefficient of friction u between the road surface and the wheel or 

45 of the braking torque or of the braking force against the wheel with respect to the slip speed. Accordingly, when these 
gradients are sleep, the peripheral speed of the wheel is great and it is hard to slip. On the contrary, when these gradi- 
ents are gentle, the peripheral speed of the wheel is small and it is easy to slip. Accordingly, these gradients express 
the slipperiness of the wheel. 

In this case, a plurality of parameters for expressing the second physical amount, for example, the relation between 

so the slip speed mentioned above and the coefficient of friction p mentioned above, has a swirl of a small hysteresis as 
the first physical amount, for example, the braking force or the like increases as mentioned above. Accordingly, the sec- 
ond physical amount detected by the second detecting means greatly varies in correspondence to the change of the 
first physical amount. 

Considering the fact that the variation of this kind is caused by the change of the first physical amount, in order to 
55 reduce the variation, it is sufficient to reduce it in correspondence to the change of the first physical amount 

Then, the reduction processing means in accordance with the present invention is structured such as to reduce the 
variation in the second physical amount detected by the second detecting means in correspondence to the change of 
the first physical amount detected by the first detecting means. In this case, the reduction process may be performed. 
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for example, by filtering at each time when the first physical amount changes by a predetermined amount or by averag- 
ing the second physical amount while the first physical amount changes by a predetermined amount 

As mentioned above, since the variation in the second physical amount expressing the slipperiness of the wheel 
due to the change of the first physical amount expressing the braking concfition or the driving condition of the wheel in 
5 the vehicle is reduced in correspondence to the change of the first physical amount, it is possible to reduce the variation 
in the second physical amount. 

In this case, since the second physical amount mentioned above is used in the road surface condition estimating 
apparatus, the ABS apparatus, the driving force controlling apparatus and the like, the third invention mentioned above 
can be also applied to the road surface condition estimating apparatus, the ABS apparatus, the driving force controlling 
w apparatus and the like. 

An invention to which the third invention estimating apparatus will be explained as follows. 

The present invention comprises first detecting means for detecting a first physical amount expressing a brake con- 
dition or a driving condition of a wheel of a vehicle; 

75 second detecting means for detecting a vehicle speed of said vehicle; 

third detecting means for detecting a second physical amount which represents slipperiness and which is deter- 
mined in accordance with a road surface condition on which said wheel is running, said first physical amount, and 
said vehicle; 

reduction processing means for reducing variations in the second physical amount detected by said third detecting 
20 means in accordance with a change in said first physical amount; 

calculating means for calculating the second physical amount for each of a plurality of predetermined road surface 
concfitions on the basis of the first physical amount and the vehicle speed; 

comparing means for comparing a detected track of the second physical amount in which variation is reduced by 
said reduction processing means with respect to the first physical amount with each of calculated tracks of the see- 
as ond physical amount calculated for each of the plurality of road surface conditions by said calculating means with 
respect to said first physical amount; and 

estimating means for estimating a road surface condition corresponding to a track closest to said detected track 
among said calculated tracks as an actual road surface condition on which said wheel is running on the basis of 
the comparison result of said comparing means. 

30 

The first detecting means detects the first physical amount expressing the braking condition or the driving condition 
of the wheel of the vehicle, the second detecting means detects the speed of the vehicle (the vehicle speed) and the 
third detecting means detects the second physical amount expressing the slipperiness of the wheel. The second phys- 
ical amount is defined by the road surface condition in which the wheel is running, the first physical amount and the 

35 vehicle speed. Concretely speaking, the second physical amount is a gradient of the coefficient of friction \i between 
the road surface and the wheel or the wheel, or a gradient of the braking torque or of the braking force with respect to 
the slip speed, as mentioned above. 

The reduction processing means reduces the variation in the second physical amount in correspondence to the 
change of the first physical amount detected by the first detecting means, as mentioned above. 

40 The calculating means calculates the second physical amount for each of a plurality of predetermined road surface 
conditions on the basis of the first physical amount detected by the first detecting means and the vehicle speed detected 
by the second detecting means. The calculating means may previously calculate the relation defined in correspondence 
to the vehicle speed between the first physical amount and the third physical amount for calculating the second physical 
amount for each of a plurality of road surface conditions, may calculate the third physical amount for each of a plurality 

45 of road surface conditions on the basis of the detected first physical amount and the detected vehicle speed, and may 
calculate the second physical amount for each of a plurality of road surface conditions from the calculated third physical 
amount Further, the calculating means may previously calculate the relation defined in correspondence to the vehicle 
speed between the first physical amount and the second physical amount for each of a plurality of road surface condi- 
tions, and may calculate the second physical amount for each of a plurality of road surface conditions on the basis of 

so the detected first physical amount and the vehicle speed. 

The comparing means compares the track of the second physical amount in which the variation is reduced by the 
reduction processing means with respect to the first physical amount expressing the braking condition or the driving 
condition of the wheel with each of the tracks of the second physical amounts calculated for each of a plurality of road 
surface conditions by the calculating means with respect to the first physical amount. The estimating means estimates 

55 the road surface condition corresponding to the track closest to the detected track among the calculated tracks as the 
actual road surface condition in which the wheel is running, on the basis of the comparison result of the comparing 
means. 

As mentioned above, the second physical amount in which the variation is reduced in correspondence to the 
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change of the first physical amount is detected and the track of the detected second physical amount is compared with 
each of the tracks of the calculated second physical amount calculated in correspondence to each of a plurality of road 
surface conditions. Accordingly, the road surface condition corresponding to the selected track corresponds to the 
actual road surface condition in which the wheel is running when selecting the track closest to the detected track among 
5 the tracks calculated in correspondence to each of a plurality of road surface conditions on the basis of the comparison 
result. Thus, it is possible to accurately estimate the actual road surface condition. 

In this case, the second physical amount mentioned above includes an error due to various causes. 

The first detecting means, the second detecting means and the third detecting means respectively detect the first 
physical amount the vehicle speed and the second physical amount at each predetermined time from the time the brak- 
10 ing is started against the wheel, and the calculating means calculates the second physical amount for each of a plurality 
of road surface conditions at each time the first physical amount changes by a predetermined amount from the time the 
braking is started against the wheel. 

Then, at each time the first physical amount changes by a predetermined amount from the time the braking is 
started against the wheel, the comparing means may compare the second physical amount detected at the predeter- 
75 mined time with each of the second physical amounts calculated for each of the plurality of road surface conditions at 
each time the first physical amount changes by a predetermined amount. 

As a result, it is possible to reduce an influence due to the error, so that it is possible to more accurately estimate 
the road surface condition. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram which shows a schematic structure of a road surface condition estimating apparatus in 
accordance with the present invention (corresponding to a first embodiment). 

Figs. 2A to 2C are graphs which show a correlation between wheel motion amounts expressing a characteristic of 
25 a road surface \i, in which Fig. 2A shows a relation between a slip ratio and the road surface fi. Fig. 2B shows a relation 
between the slip ratio and a braking force and Fig. 2C shows a relation between the braking force for each of the vehicle 
speeds and a gradient Qd of the road surface n- 

Fig. 3 is a block diagram which shows a schematic structure of a road surface condition estimating apparatus in 
accordance with the present invention (corresponding to a second embodiment). 
30 Fig. 4 is a graph which shews a relation between a braking force P c and a slip speed Av for each of the road surface 
conditions and for each of the vehicle speeds. 

Figs. 5A and 5B are graphs which shew a relation between a slip speed Av and a gradient of the road surface 
\i for each road surface condition and for each vehicle speed. Fig. 5A is a graph for a dry road surface and Fig. 5B is a 
graph for a snow covered road surface. 
35 Figs. 6A and 6B are graphs which show a relation between a braking force P c and the gradient Gy of the road sur- 
face (i in the case of setting the slip speed Av as a parameter for each of the vehicle speeds. Fig. 6 A is a graph for a dry 
road surface and Fig. 6B is a graph for a snow covered road surface. 

Figs. 7A to 7C are graphs which show a relation between the braking force P c and the gradient of the road sur- 
face n for each vehicle speed and road surface condition obtained on the basis of Figs. 6A and 6B. Fig. 7A is a graph 
40 in the case of the vehicle speed V = 20 km/h, Fig. 7B is a graph in the case of the vehicle speed V = 40 km/h, and Fig. 
7C is a graph in the case of the vehicle speed V = 60 km/h. 

Fig. 8 is a block diagram which shows a schematic structure of a road surface condition estimating apparatus in 
accordance with the present invention (corresponding to a third embodiment). 

Figs. 9A and 9B are graphs which show the maximum braking force for each vehicle speed and road surface con- 
45 ditions. Fig. 9A is a graph in the case of the vehicle speed V = 40 km/h and Fig. 9B is a graph in the case of the vehicle 
speed V = 60 km/h. 

Rg. 10 is a block diagram which shows a structure of the road surface condition estimating apparatus in accord- 
ance with the first embodiment of the present invention (a first example). 

Fig. 11 is a Wock diagram which shows a structure of the road surface condition estimating apparatus in accord- 
50 ance with the first embodiment of the present invention (a second example). 

Rg. 12 is a block diagram which shows a structure of the road surface condition estimating apparatus in accord- 
ance with the second embodiment of the present invention. 

Rg. 13 is a block diagram which shows a structure of the road surface condition estimating apparatus in accord- 
ance with the third embodiment of the present invention. 
55 Rg. 14 is a block diagram which shows a structure of a peak n searching portion in accordance with an embodi- 
ment of the present invention. 

Rg. 15 is a block diagram which shows a structure of a resonance gain calculating portion in accordance with an 
embodiment of the present invention. 
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Fig. 1 6 is a schematic view which shows an equivalent model of a vibration system constituted by a wheel, a vehicle 
body and a road surface. 

Fig. 17 is a schematic view which shews a change characteristic of a change in a coefficient of friction u with 
respect to the slip speed and in order to explain that a resonance gain is equivalent to a braking torque gradient, also 
5 shows that a change of n around a center of a micro vibration can be approximated by a straight line. 
Fig. 18 is a block diagram which shows a hardware structure of a brake portion. 

Fig. 19 is a schematic view which shows a command to a control valve in the case of simultaneously performing 
micro vibration of a braking pressure and control of an average braking force. 

Fig. 20 is a schematic view which shows an embodiment of a system construction in the case of calculating a 
10 parameter showing a friction condition from a calculated road surface condition and a slip ratio in the second embodi- 
ment of the present invention. 

Fig. 21 is a graph which shows a characteristic of a coefficient of friction ji between a tire and a road surface with 
respect to a slip ratio X. 

Fig. 22 is a block diagram which shows a road surface condition estimating apparatus in accordance with a fourth 
is embodiment. 

Fig. 23 is a block diagram which shows a wheel speed micro amplitude detecting portion. 

Fig. 24 is a graph which shows a relation between a braking force and the slip speed at a predetermined speed in 
a Dry road surface condition. 

Fig. 25 is a graph which shows a relation between a braking force and the slip speed at a predetermined speed in 
20 a Snow covered road surface condition. 

Fig. 26 is a graph which shows a relation between a braking force and the slip speed at a predetermined speed in 
an Ice covered road surface condition. 

Fig. 27 is a graph which shows a relation between a braking force and the road surface nat a predetermined speed 
in a Dry road surface condition. 
25 Fig. 28 is a graph which shows a relation between a braking force and the road surface fiat a predetermined speed 
in a Snow covered road surface condition. 

Fig. 29 is a graph which shows a relation between a braking force and the road surface fiat a predetermined speed 
in an Ice covered road surface condition. 

Fig. 30 is a flow chart which shows an identification processing routine which a time series data track identifying 
30 portion performs at a time of starting the brake. 

Fig. 31 is a schematic view which shows a state of identification process (track) of the time series data track iden- 
tifying portion. 

Fig. 32 is a graph which shows a relation between the braking force and the gracfient of the road surface \x. 
Fig. 33 is a graph which shews a relation between a braking force and the gradient the road surface u at a prede- 
35 termined speed in a Dry road surface condition. 

Fig. 34 is a graph which shows a relation between a braking force and the gradient of the road surface ji at a pre- 
determined speed in a Snow covered road surface condition. 

Fig. 35 is a graph which shows a relation between a braking force and the gradient of the road surface fiat a pre- 
determined speed in an Ice covered road surface condition. 
40 Fig. 36 is a block diagram which shows a road surface condition estimating apparatus in accordance with a fifth 
embodiment. 

Fig. 37 is a graph which shows a relation between a braking force and the gradient of the road surface n at a pre- 
determined speed in a Dry road surface condition. 

Fig. 38 is a graph which shows a relation between a braking force and the gradient of the road surface \i at a pre- 
45 determined speed in a Snow covered road surface condition. 

Figs. 39A and 39B are graphs which respectively show characteristic of a coefficient of friction \i between a tire and 
a road surface and that of the gradient of the road surface ji with respect to a slip speed Av. 

Fig. 40 is a graph which shows a characteristic of the coefficient of friction \i between the tire and the road surface 
with respect to the slip speed in a micro manner. 
so Fig. 41 is a graph which shows a characteristic of the resonance gain calculated in the resonance gain calculating 
portion with respect to the braking force. 

Fig. 42 is a graph which shows a resuft obtained by filtering the resonance gain calculated in the resonance gain 
calculating portion for each fixed sampling times as a characteristic to the braking force. 

Fig. 43 is a graph which shows a result obtained by filtering the resonance gain calculated in the resonance gain 
55 calculating portion for each time when the braking force changes at a predetermined amount as a characteristic to the 
braking force, and compared with the result in Fig. 42. 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Each of embodiments of a road surface condition estimating apparatus in accordance with the present invention 
will be described in detail below with reference to the drawings. Here, in the embodiments, it is supposed that the road 
5 surface condition estimating apparatus is applied to a vehicle. 

(First Embodiment) 

Fig. 1 0 shows a block diagram of a structure of the road surface condition estimating apparatus in accordance with 
10 the first embodiment of the present invention. In this case, the structure of the first embodiment is obtained by embod- 
ying the schematic structure shown in Fig. 1 . 

As shown in Fig. 10, the road surface condition estimating apparatus 10 in accordance with the first embodiment 
is provided with converting portions 12, 14 and 16 which respectively store a dji/dAv table showing a relation among 
three wheel motion amounts, that is, a gradient d\JdAv (which is supposed to be a resonance gain in the case of this 
is embodiment) of a road surface *i with respect to a slip speed Av (or a slip ratio), a braking force P c , and a vehicle speed 
V for each of the road surface conditions, i.e. a dry road surface (Dry), a snow covered road surface (Snow) and an ice 
covered road surface (Ice). 

When the vehicle speed and the braking force are input, the converting portions 12, 14 and 16 respectively output 
resonance gains Gy 1( and Q43 (although in Fig. 10, a mark A showing an estimation is attached, in the following 
20 description the mark A is omitted) corresponding to the input values on the basis of a table (refer to Fig. 2C) in connec- 
tion to a function stored therein. 

In this case, the converting portions can be constructd as a ROM which previously stores all the converted values 
with respect to the input values in a predetermined range, and outputs the converted value corresponding to the data 
value input thereto. Further, the converting portions can be constituted as a neutral network which can learn a non-lin- 
25 ear function such as the d\ifdAv table, an application filter or the like. 

Further, the road surface condition estimating apparatus 10 is provided with subtracting devices 18, 20 and 22 for 
respectively calculating differences - Qdi. Qd ' and Qd ' by respectively subtracting a resonance gain G^, 
Gd2 and converted by the converting portions 12, 14 and 16 from a detected value Qd of a resonance gain input to 
the apparatus 10. 

30 Still further, the road surface condition estimating apparatus 10 is provided with a minimum value selecting portion 
24 for selecting a minimum value among absolute vales of the differences respectively calculated by the subtracting 
devices 18, 20 and 22. The minimum value selecting portion 24 specifies a converting portion outputting the resonance 
gain which gives the minimum value, and outputs an information of the road surface condition corresponding to this 
converting portion. 

35 Furthermore, a vehicle speed detecting portion 26 for detecting vehicle speed, a braking force detecting portion 34 
for detecting a braking force acting on the wheel and a resonance gain calculating portion 36 as road surface \i gradient 
detecting means in Fig. 1 are provided in an outer portion of the road surface condition estimating apparatus 10. The 
vehicle speed detecting portion 26 and the braking force detecting portion 34 are respectively connected to the convert- 
ing portions 12, 14 and 16, and the resonance gain calculating portion 36 is respectively connected to the subtracting 

40 devices 18, 20 and 22. 

Among them, the vehicle speed delecting portion 26 is constituted by a wheel speed sensor 28 for detecting the 
wheel speed, a brake start judging portion 30 for judging a time when the driver presses the brake pedal as the brake 
start point, and a brake start vehicle speed determining portion 32 for outputting a wheel speed V w detected at the 
judged brake start point as vehicle speed V. Thus, in the vehicle speed detecting portion 26 in accordance with this 

45 embodiment, the vehicle speed is estimated by supposing that the wheel speed at a time of starting the brake substan- 
tially equals to the vehicle speed. Alternatively, the vehicle speed detecting portion 26 may be structured as a vehicle 
speed sensor for directly detecting the vehicle speed. 

Further, the braking force detecting portion 34 estimates the braking force acting on the wheel from the road sur- 
face as a friction force in the following manner in accordance with a mechanical model of the wheel. 

so A braking torque T B acting on the wheel to an opposite direction to the rotating direction of the wheel and a tire 
torque T f applied by a braking force F acting on the wheel as a friction force to a rotating direction of the wheel act on 
the wheel. The braking torque T B is caused by a braking force acting on the brake disc of the wheel in such a manner 
as to prevent the wheel from rotating, and the braking force F and the tire torque T f can be expressed by the following 
formula when a coefficient of friction between the wheel and the road surface is set to be ^ a radius of the wheel is 

55 set to be rand a wheel load is set to be W. 

F = ji B W 
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T f = F x r= u B Wr 



Accordingly, an equation of motion of the wheel is expressed as follows. 



l^ = H B Wr-T B = Fr-T B (15) 



10 in which reference symbol I is a moment of inertia of the wheel and © is a rotational speed of the wheel(a wheel 

speed). 

By detecting the wheel acceleration (d©/dt) and calculating the braking torque T B on the basis of a wheel cylinder 
pressure applied to the brake disc, the braking force F can be estimated on the basis of the formula (15). Concretely 
speaking, it is constructed by setting a driving torque of the wheel calculated from an acceleration opening degree and 
is the like, and a mechanical model equivalent to the formula (15) in which the braking force F as a disturbance acts on 
the wheel to be an observer. In this observer, the disturbance of the equivalent model and the rotational speed are mod- 
ified at every control periods so that a deviation between the rotational position obtained by integrating the formula (15) 
for two stages and the rotational position actually detected equal zero, and the modified disturbance is estimated as the 
braking force. 

20 Further, the resonance gain calculating portion 36 is constituted by a wheel speed micro amplitude detecting por- 
tion 40 for detecting a micro amplitude of a resonance frequency ©oo of the wheel speed V w (a wheel speed micro 
amplitude ©wv) at a time of minutely vibrating the braking pressure at a resonance frequency <doo (the formula (3)) in 
the vibration system comprising the vehicle body, the wheel and the road surface around an average braking pressure, 
braking pressure micro amplitude detecting portion 42 for detecting a micro amplitude P v of the braking pressure at the 

25 resonance frequency ©qo, and a dividing device 44 for outputting a resonance gain obtained by dividing a detected 
wheel speed micro amplitude ©^ by the braking pressure micro amplitude P v as shown in Fig. 15. 

Among them, the wheel speed micro amplitude detecting portion 40 can be realized as a calculating portion for per- 
forming a filter treatment for extracting a vibration component of the resonance frequency <dqo. For example, since the 
resonance frequency ©oo of the vibration system is about 40 Hz, a period is set to be 24 ms and about 41 .7 Hz consid- 

30 ering controllability the above, a band-pass filter having a central frequency of this frequency is provided, and full wave 
of the filter output are rectified and made smooth to a direct current, thereby outputting the wheel speed micro ampli- 
tude. Further, the structure can be realized so as to continuously take in the time series data of integral times of the 
period, for example, 24 ms of a period and 48 ms of two periods and calculate the correlation between a unit sine wave 
and a unit cosine wave of 41 .7 Hz. 

35 Here, micro vibration means for applying the braking pressure micro amplitude P v having the resonance frequency 
around an average braking pressure P m will be described below. At first, a portion for converting the average braking 
pressure command and the micro vibration command into the actual braking torque to a wheel (a valve control system) 
comprises a master cylinder 48, a control valve 52, a wheel cylinder 56, a reservoir 58 and an oil pump 60, as shown 
in Fig. 18. 

40 Among them, the brake pedal 46 is connected to a pressure increasing valve 50 of the control valve 52 via the mas- 
ter cylinder 48 which increases a pressure in correspondence to a stepping force applied to the brake pedal 46. Further, 
the control valve 52 is connected to the reservoir 58 serving as a tow pressure source via the pressure reducing valve 
54. Still further, a wheel cylinder 56 for applying a braking pressure supplied by the control valve to the brake disc is 
connected to the control valve 52. The control valve 52 controls opening and closing of the pressure increasing valve 

45 50 and the pressure reducing valve 54 on the basis of an input valve operation command. 

In this case, when the control valve 52 is controlled to open only the pressure increasing valve 50, a hydraulic pres- 
sure of the wheel cylinder 56 (a wheel cylinder pressure) is increased to a hydraulic pressure of the master cylinder 48 
(a master cylinder pressure) which is in proportion to the pressure obtained by stepping the brake pedal 46 by the 
driver. On the contrary, when the control is performed so as to open only the pressure reducing valve 54, the wheel cyl- 

50 inder pressure is decreased to a pressure of the reservoir 58 (a reservoir pressure) which substantially corresponds to 
an atmospheric pressure. Further, when the control is made so as to close both of the valves, the wheel cylinder pres- 
sure is maintained. 

The braking force applied to the brake disc from the wheel cylinder 56 (corresponding to the wheel cylinder pres- 
sure) can be calculated by a rate among a pressure increasing time for supplying a high hydraulic pressure of the mas- 
55 ter cylinder 48, a pressure reducing time for supplying a low hydraulic pressure of the reservoir 58 and a maintaining 
time for holding the supply hydraulic pressure, the master cylinder pressure detected by the pressure sensor or the like, 
and the reservoir pressure. 

Accordingly, a desired braking torque can be realized by controlling the pressure increasing and reducing times of 
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the control valve 52 in accordance with the master cylinder pressure. Then, a micro vibration of the braking pressure 
can be realized by performing the pressure increasing and reducing control at a period corresponding to the resonance 
frequency at the same time of the pressure increasing and reducing control of the control valve 52 realizing the average 
braking force. 

5 As the contents of a concrete control, as shown in Fig. 19, respective modes comprising a pressure increasing 
mode and a pressure reducing mode are switched at every half period T/2 of the period of the micro vibration (for exam- 
ple, 24 ms), a pressure increasing command and a pressure reducing command are output to the valve for the respec- 
tive time of the pressure increasing time tj and the pressure reducing time t r from a moment of switching the mode, and 
a pressure holding command is output for the remaining time. The average braking force is defined by the rate between 

10 the pressure increasing time tj and the pressure reducing time ^ in correspondence to the master cylinder pressure, and 
the micro vibration is applied around the average braking force by switching the pressure increasing and pressure 
reducing modes at every half period T/2 corresponding to the resonance frequency. 

In this case, since the braking pressure micro amplitude P v is defined by the master cylinder pressure, a time length 
of the pressure increasing time tj and the time length of the pressure reducing time ^ of the valve shown in Fig. 19 in 

75 accordance with a predetermined relation, the braking pressure micro amplitude detecting portion 42 shown in Fig. 1 5 
can be structured as a table for outpirtting the braking pressure micro amplitude P v from the master cylinder pressure, 
the pressure increasing time tj and the pressure reducing time t r 

Further, the road surface condition estimating apparatus 10 in accordance with the first embodiment is provided 
with a peak n searching portion 61 as shown in Fig. 14. The peak y. searching portion 61 is constituted by a memory 

20 64 for storing the data of the peak \i value for every road surface condition and the slip ratio giving the peak u. value, 
and a data searching portion 62 for searching the data in the memory 64 on the basis of an address table 63 showing 
an address of the data for each of the road surface conditions. 

Next an operation of the first embodiment will be described below. 

When the vehicle speed V detected by the vehicle speed detecting portion 26 and the braking force P c detected by 
25 the braking force detecting portion 34 are input to each of the converting portions 12,14 and 1 6 of the road surface con- 
dition estimating apparatus 10, each of the converting portions converts two input detected values in accordance with 
the d\ddAv table given to each of them, and outputs each of the resonance gains Qdi. and Gd3 corresponding to 
these detected values. 

When the resonance gain calculated by the resonance gain calculating portion 36 is respectively input to the 
30 subtracting devices 18, 1 9 and 20 together with the respective converted resonance gains , and G^, these sub- 
tracting devices respectively calculate the differences Gd - Gdi. Gd -G<& and Gd - Gd3 and output them to the minimum 
value selecting portion 24. 

The minimum value selecting portion 24 selects the difference having the smallest absolute value among the cal- 
culated differences, and specifies the converting portion which gives the selected minimum absolute difference. Then, 

as the road surface condition corresponding to the specified converting portion is output as an information code of a pre- 
determined type or the like. The information code expressed by the type corresponds to any of Dry, Snow and Ice, and 
is expressed by a type which the control portion of the vehicle employing the road surface condition estimating appara- 
tus 10 can recognize. For example, in the case that Qd - Gd3 becomes the minimum value, the information code having 
the road surface condition of "Ice" is output. 

40 Since the d^i/dAv table expressing the correlation among the vehicle speed, the braking force and the resonance 
gain precisely shows the road surface \i characteristic, the road surface condition can be estimated in a significantly 
accurate manner by comparing the relation between the respective wheel motion amounts which the du/dAv table 
shews in accordance with the present invention with the relation between the actual detected values. Further, in accord- 
ance with this embodiment, since the resonance gain Gd sensitively expressing the friction state on the basis of the res- 

45 onance characteristic in the vibration system is employed, in comparison with the case of using other wheel motion 
amounts, the road surface condition can be more accurately estimated. 

Further, in the peak \i searching portion 61 , the data searching portion 62 interprets the information code of the 
road surface condition output by the minimum value selecting portion 24 and reads out the address corresponding to 
the information code from the address table 63. Then, the data searching portion 62 searches the peak \i value and the 

so slip ratio of the corresponding address in the memory 64 and outputs them. Accordingly, the output data becomes the 
peak ii value of the road surface on which the vehicle currently runs, and the slip ratio which gives the peak \i value on 
the road surface. 

As mentioned above, in accordance with this embodiment, since the road surface condition of the road surface on 
which the vehicle currently runs and the peak \i value can be accurately estimated, safe driving can be realized regard- 
55 less of the road surface condition by employing the estimated result to the vehicle stabilizing control such as VSC, ABS, 
TRC and the like and the alarm to the driver. 

Here, in the embodiment shown in Rg. 10, the vehicle speed, the braking force and the resonance gain are used 
as the wheel motion amounts. The present invention is not limited to this, however, for example, the structure may be 
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made in a manner as shown in Fig. 1 1 . 

In Fig. 11 , the slip speed is used as one of the wheel motion amounts in place of the braking force. Accordingly, 

each of the converting portions 13,15 and 1 7 shown in Fig. 1 1 is provided with a dji/dAv table showing a relation among 

the vehicle speed, the slip speed and the resonance gain. Then, in place of the braking force detecting portion 34, there 
5 is provided a slip speed calculating device 38 for calculating the slip speed from the wheel speed detected by the wheel 

speed sensor 28 and the vehicle speed determined by the brake start vehicle speed determining portion 32, and the 

slip speed calculating device 38 is connected to each of the converting portions 13, 15 and 17. 

Also in the structure shown in Fig. 11, since the road surface \i characteristic can be expressed by the slip speed, 

the vehicle speed and the resonance gain, the road surface condition can be accurately estimated. Further, since the 
10 braking force detecting portion 34 is omitted and it is sufficient to provide the simple slip speed calculating device in 

place thereof, the whole apparatus can be made simple. 

(Second Embodiment) 

is Fig. 1 2 shows a block diagram of a structure of a road surface condition estimating apparatus in accordance with 
a second embodiment of the present invention embodying the schematic structure shown in Fig. 3. In this case, the 
same reference numerals are attached to the same elements as those in the first embodiment and the detailed descrip- 
tion thereof will be omitted. 

As shown in Fig. 1 2, a road surface condition estimating apparatus 10b in accordance with the second embocfiment 
20 is constituted by a vehicle speed calculating portion 70 for calculating the vehicle speed V from the wheel speed V w 
detected by the wheel speed sensor 65 and the braking force P c detected by the braking force detecting portion 34, and 
a road surface condition calculating portion 80 for calculating the road surface condition from the vehicle speed V cal- 
culated by the vehicle speed calculating portion 70, the detected braking force P c and the resonance gain Qd. 

Among them, the vehicle speed calculating portion 70 is provided with converting portions 71 , 72, 73, ... for respec- 
ts tively storing a Av table showing a relation between the braking force and the resonance gain at each of fixed sip 
speeds (0.1 m/s, 0.4 m/s, 0.8 m/s, ... ). 

The converting portions 71, 72, 73, ... respectively output the resonance gains G^-j, G^, G^, ... (although in Fig. 
1 2, a mark A showing an estimation is attached, the mark A is omitted in the following description) corresponding to the 
input value on the basis of the Av table (refer to Figs. 6A and 6B) with respect to the function stored therein when the 
30 detected braking force is input 

Further, the vehicle speed estimating portion 70 is provided with subtracting devices 74, 75 and 76 for respectively 
calculating differences Qd - Qdi» Qd - Qd2» Qd " Qd3> - respectively subtracting a resonance gain Gdi, Q^, Q^, ... 
converted by the converting portions 71 , 72, 73, ... from a detected value Qd of a resonance gain input to the apparatus 
106. 

35 Still further, the vehicle speed estimating portion, 70 is provided with a minimum value selecting portion 77 for 
selecting a minimum value among absolute vales of the differences respective calculated by the subtracting devices 74, 
75, 76, ... . The minimum value selecting portion 77 specifies a converting portion outputting the resonance gain which 
gives the minimum value, and outputs a slip speed Av corresponding to the converting portion. 

Furthermore, a vehicle speed calculating device 78 is connected to the minimum value selecting portion 77. The 

40 vehicle speed calculating device 78 calculates the vehicle speed V by adding the wheel speed V w detected by the wheel 
speed sensor 65 to the slip speed Av output from the minimum value selecting portion 77. 

In this case, as shown in Fig. 12, it is also possible to output the slip speed Av output from the minimum value 
selecting portion 77 and the vehicle speed V calculated by the vehicle speed calculating device 78 to an outer portion 
of the road surface condition estimating apparatus 10b. Further, it is also possible to provide a dividing device 81 for 

45 outputting the slip ratio X by dividing the slip speed Av output from the minimum value selecting portion 77 by the vehicle 
speed V calculated by the vehicle speed calculating device 78. 

The road surface condition calculating portion 81 in accordance with the second embodiment is structured in the 
same manner as the road surface condition estimating apparatus 10 shown in Fig 10 except for using the vehicle speed 
calculated by the vehicle speed calculating portion 70, not by the vehicle speed detecting portion 26. Accordingly, the 

so detailed description thereof will be omitted. In this case, it is possible to structure the road surface condition calculating 
portion 81 in such a manner as to use the vehicle speed and the slip speed as in the road surface condition estimating 
apparatus 1 0 shown in Fig. 1 1 and to input the vehicle speed V and the slip speed Av calculated by the vehicle speed 
calculating portion 70 into the road surface condition calculating portion 81. 
Next an operation of the second embodiment will be described below. 

55 When the braking force P c detected by the braking force detecting portion 34 is input to each of the converting por- 
tions 71, 72, 73, ... of the vehicle speed calculating portion 70, each of the converting portions converts the input 
detected values in accordance with the Av table given to each of them, and outputs each of the resonance gains Q d1 , 
Qd2. Qd3. ... corresponding to the detected values. 
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When the resonance gain calculated by the resonance gain calculating portion 36 are respectively input to the 
subtracting devices 74, 75, 76, ... together with the converted resonance gains G^, G^, .... these subtracting 
devices respectively calculate the differences - Gd 1( - G^, Gd ... and output them to the minimum value 
selecting portion 77. 

5 The minimum value selecting portion 77 selects the difference having the smallest absolute value among the cal- 
culated differences, and specifies the converting portion which gives the selected minimum absolute difference. Then, 
the slip speed Av correspond ng to the specified converting portion is output. For example, in the case that Gd - 
becomes the minimum value, the slip speed is output as 0.4 m/s corresponding to the converting portion 72. 

Then, the vehicle speed calculating device 78 calculates the vehicle speed V by adding the slip speed Av output by 

w the minimum value selecting portion 77 and the wheel speed V w detected by the wheel speed sensor 65, and outputs 
it to the road surface condition calculating portion 80. In the road surface condition calculating portion 80, when the 
vehicle speed V and the braking force is input the road surface concGtion of the road surface on which the vehicle cur- 
rently runs is calculated and output in the same manner as that of the first embodiment. 

Further, the slip speed Av output by the minimum value selecting portion 77 and the vehicle speed V calculated by 

75 the vehicle speed calculating device 78 are output to the outer portion of the apparatus, and further Av/V is calculated 
by the dividing device 81 , so that the calculated result is output as the slip ratio X. 

As mentioned above, in accordance with the second embodiment, in spite of having the dp/dAv table including the 
vehicle speed as one of the wheel motion amounts, it is not necessary to directly detect the vehicle speed, and the 
wheel speed is detected in place thereof. Since it is hard to accurately estimate the vehicle speed in comparison with 

20 the wheel speed, it is possible to more accurately estimate the road surface condition by estimating the vehicle speed 
on the basis of the Av table and the wheel speed as in this embodiment in comparison with the first embocfiment using 
the detected vehicle speed. 

Particularly, in the embodiment having the vehicle speed detecting portion 26 shown in Fig. 10, since it is supposed 
that the wheel speed at a time of starting brake substantially equals to the vehicle speed, the road surface condition can 
25 be accurately calculated only as far as the supposition is established. On the other hand, in this embodiment it is pos- 
sible to obtain the vehicle speed as accurately as possible by providing a detailed Av table. 

Further, in this embodiment, since the accurate slip ratio is obtained by accurately estimating the vehicle speed, 
various application can be performed by using the slip ratio. For example, as shown in Ftg. 20, the peak ji searching 
portion 61 searches the peak pofa'the corresponding road surface and the slip ratio Xq giving the peak no by using the 
30 road surface condition calculated by the road surface condition estimating apparatus 10b. 

In the case of applying this embodiment to the anti-lock brake control, the braking force is controlled so that the cur- 
rent slip ratio X calculated by the road surface condition estimating apparatus 10b coincides with the searched slip ratio 
Xq. Accordingly, since the calculated value of the slip ratio Ao is changed in correspondence to the change of the road 
surface condition, a stable anti-lock brake control can be always obtained without relation to the road surface condition. 
35 Further, it is possible to provide a road surface \i calculating portion 82 having the road surface ji table showing the 
relation of the road surface ji with respect to the slip ratio for each of the road surface conditions, which selects the road 
surface \i table corresponding to the input road surface condition and outputs the p value corresponding to the calcu- 
lated slip ratio in accordance with the selected table, in the system shown in Ftg. 20. 

When the road surface \i value of the road surface on which the vehicle currently runs is obtained by the road sur- 
40 face \i calculating portion 82, it is possible to estimate the area in the road surface \i characteristic to which the current 
friction state belongs by comparing the road surface \i value with the peak no value in the corresponding road surface. 

As a parameter showing the friction state, for example, jiq ~ M can 00 used The wheel may be undesirably locked 
when the parameter frio - n) becomes a negative value. Therefore, a system for giving an alarm to the driver when it 
happens, results in safe driving. In this case, Xq - X may be employed as a parameter for expressing the friction state of 
45 this kind. 

In accordance with this embodiment, since not only the road surface condition but also the peak \i value and the 
current friction condition can be known, a limit value of the braking force (driving force) can be clarified. Accordingly, in 
the vehicle stabilizing control, more stable control can be realized by changing the control gain and setting the control 
target value on the basis of these parameters. 

50 

(Third Embodiment) 

Fig. 13 shows a block diagram of a structure of a road surface condition estimating apparatus in accordance with 
a third embodiment of the present invention which embodies the schematic structure of Fig. 8. In this case, the same 
55 reference numerals are attached to the same elements as those of the first and second embodiments, and the detailed 
description thereof will be omitted. 

As shown in Fig. 13, a road surface condition estimating apparatus 10c in accordance with the third embodiment 
is constituted by a road surface condition estimating portion 83 for estimating the road surface condition and the vehicle 
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speed V in accordance with the detected wheel speed the resonance gain and the braking force P c , a maximum 
braking force selecting portion 84 for calculating a maximum braking force against the road surface on which the vehicle 
currently runs in accordance with the estimated road surface condition and the vehicle speed, and a subtracting device 
89 for outputting a surplus degree AP C of the braking force by subtracting the detected braking force from the calculated 
5 maximum braking force. In this case, the surplus degree AP C of the braking force is calculated for each of the wheels. 

Among them, the maximum braking force selecting portion 84 is provided with converting portions 85, 86 and 87 
having the maximum braking force table (refer to Figs. 9A and 9B) giving the maximum braking force at each of the vehi- 
cle speeds, in each of the road surface conditions. These converting portions respectively calculate the maximum brak- 
ing force corresponding to the input vehicle speed V on the basis of the maximum braking force table, and output each 
w of the calculated maximum braking forces to a table selecting portion 88 in the following step. The table selecting por- 
tion 88 selects only the maximum braking force converted by the converting portion corresponding to the estimated 
road surface condition and outputs it. 

Further, the road surface condition estimating portion 83 can be constructed in the same manner as the road sur- 
face condition estimating apparatus 10b shown in Fig. 12. 
15 Next an operation of the third embodiment will be described below. 

In the third embodiment, since the surplus degree AP C of the braking force is calculated for each of the wheels, a 
spin or the like is prevented, for example, by giving larger braking force and driving force to the wheel having a high sur- 
plus degree and controlling the braking force of each of the wheel to be equal to each other, so that a more stable vehi- 
cle running can be realized. 

20 Although three of the embodiments in accordance with the present invention are described above, the present 
invention is not limited to the embodiments mentioned above, and can be optionally and preferably modified within the 
range of the present invention. 

For example, in the embodiment, the road surface condition is classified into three kinds comprising a dry road sur- 
face, a snow covered road surface and an ice covered road surface. The other kinds of road surfaces, however, such 

25 as a gravel road, a sand road surface, a wet road surface and the like can be estimated, and the road surface conditions 
more than three kinds can be more accurately estimated. 

Further, in Figs. 10 and 1 1 , the vehicle speed and the braking force are converted into the resonance gain using 
each of the tables and the converted value of the resonance gain is compared with the detected value thereof. Other 
than that, the structure may be made such that the resonance gain and the vehicle speed, or the resonance gain and 

30 the braking force are converted into the braking force or the vehicle speed in accordance with each of the tables, and 
the converted value and the detected value of the same wheel motion amount are compared with each other, thereby 
estimating the road surface condition. Further, the same manner can be applied to the Av table in Fig. 12. 

Still further, in the embodiments, the resonance gain is used for expressing the gradient of the road surface \l Alter- 
natively, the braking torque gradient may be employed. 

35 Next, a fourth embodiment in accordance with the present invention will be described in detail with reference to the 
drawings. 

As shown in Fig. 22, a road surface condition estimating apparatus in accordance with this embodiment is provided 
with a wheel speed sensor 28 for detecting a wheel speed, a braking force detecting portion 34 for detecting a braking 
force acting on the wheel, a resonance gain calculating portion 36 for calculating a resonance gain mentioned below as 
40 a road surface \i gradient (an gradient of a coefficient of friction between the road surface and the wheel with respect to 
the slip speed), and a calculating circuit 1 0A for calculating the road surface condition and the slip speed. In this case, 
since the wheel speed sensor 28, the braking force detecting portion 34 and the resonance gain calculating portion 36 
are respectively structured in the same manner as those of the embodiments mentioned above, the detailed description 
will be omitted. 

45 The calculating circuit 1 0A is provided with reference value forming and storing circuits 1 2D, 1 2S and 1 21 for form- 
ing and storing a reference value of the road surface u. gradient corresponding to a predetermined road surface condi- 
tion, that is, Dry, Snow and Ice at every time when the braking force detected by the braking force detecting portion 34 
changes by a predetermined amount, and a resonance gain storing circuit 12J for storing the resonance gain at every 
time when the braking force changes by a predetermined amount The reference value forming and storing circuits 1 2D, 

so 1 2S and 121 and the resonance gain storing circuit 1 2J are connected to a time series data track identifying portion 1 4A 
for comparing a time series data of a reference value of the road surface \i gradient stored in the reference value form- 
ing and storing circuits 1 2D, 1 2S and 1 21 with a time series data of a reference value of the road surface \i gradient actu- 
ally detected and stored in the resonance gain storing circuit 12J and calculating the road surface condition. A slip 
speed selecting portion 16A is connected to the output end of the time series data track identifying portion 14 A. 

55 Since each of the reference value forming and storing circuits 12D, 12S and 121 is structured in the same manner, 
only the reference value forming and storing circuit 12D is explained, and the other explanations will be omitted. The 
reference value forming and storing circuit 12D is provided with a vehicle speed calculating portion 18A connected to 
the wheel speed sensor 28. A slip speed calculating portion 20A for calculating the slip speed is connected to the vehi- 
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cle speed calculating portion 18 A. The braking force detecting portion 34 is connected to the slip speed calculating por- 
tion 20A. The slip speed calculating portion 20A is connected to a ^ gradient reference value forming portion 22A for 
forming a reference value of the road surface \i gradient. A vehicle speed calculating portion 18A and the braking force 
detecting portion 34 are connected to the u gradient reference value forming portion 22A. The \i gradient reference 
5 value forming portion 22A is connected to a reference value time series data storing portion 24A for storing the time 
series data of the reference value of the road surface n gradient formed by the u gradient reference value forming por- 
tion 22 A. The reference value time series data storing portion 24A is connected to the time series data track identifying 
portion 14A. 

The resonance gain storing circuit 12J is provided with a resonance gain time series data storing portion 25 for stor- 
10 ing the time series data of the detected resonance gain. The braking force detecting portion 34 and the resonance gain 
calculating portion 36 are connected to the resonance gain time series data storing portion 25. The resonance gain 
time series data storing portion 25 is connected to the time series data track identifying portion 14A. 

Each of the slip speed calculating portions 20A of the reference value forming and storing circuits 1 2D, 1 2S and 121 
is connected to the slip speed selecting portion 16 A. 
is In this case, the wheel speed micro amplitude detecting portion 40 of the resonance gain calculating portion 36 can 
be realized as a calculating portion for performing a filter-treatment for extracting a vibration component of a resonance 
frequency of ©oo as shown in Fig. 23. For example, since the resonance frequency ©°o of this vibration system is about 
40 Hz, a period is set to be 24 ms and about 41 .7 Hz considering controllability, and a band pass filter 75 having the 
frequency as a central frequency is provided. Only the frequency component adjacent to about 41.7 Hz is extracted 
20 from a wheel speed signal <Oj by means of this filter. Further, the full wave of the filter output is rectified and smoother ed 
to a direct current by means of a full wave rectifier 76, and a wheel speed micro amplitude is output from the direct 
current smoothening signal by passing only a low band vforation component by means of a low pass filter 77. 

In this case, the wheel speed micro amplitude detecting portion 40 can be realized so as to continuously lake in the 
time series data of integral times of the period, for example, of 24 ms of a period and 48 ms of two periods and calculate 
25 the correlation between the unit sine wave and the unit cosine wave of 41 .7 Hz. 
Next, an operation of this embodiment will be explained below. 

At first a method of forming a reference value of the road surface u gradient by means of the reference value form- 
ing and storing circuits 12D, 12S and 121 will be explained. 

The gradient of the road surface \i corresponds to an gradient of the coefficient of friction \i between the road 
30 surface and the wheel with respect to the slip speed Av, and is expressed by the following formula (1 6) at each of the 
braking force, the road surface condition and the vehicle speed. 

G ^dtv^I7^K7\ (16) 

35 

in which the numbers 1 and 2 attached to the letters u and Av respectively represent the time at which the brak- 
ing changes by a predetermined amount from the time of starting the brake, wherein the number 2 is the present time, 
he number 1 is the time prior to the present time, and the number 2 will be replaced by the number 1 when the braking 
40 force changes by the predetermined amount from the present time. 

Accordingly, the gradient Gy of the road surface ^ can be specified by the braking force the road surface condi- 
tion, the vehicle speed V, the slip speed Av and the coefficient of friction jl 

In this case, it is known that the slip speed Av and the braking force P have an inherent relation in accordance with 
each of the road surface conditions and the vehicle speeds as shown in Figs. 24 to 26. 
45 Accordingly, each of the slip speed calculating portions 20A of the reference value forming and storing circuits 1 2D, 
12S and 121 stores the relational formula between the braking force P and the slip speed Av shown in Figs. 24 to 26 in 
accordance with the vehicle speed. Accordingly, when the braking force P c and the vehicle speed V are input the slip 
speed calculating portion 20A calculates the slip speed Av corresponding to the braking force P c and the vehicle speed 
V on the basis of the relational formula mentioned above, the braking force P c and the vehicle speed V, and outputs the 
so calculated slip speed Av to the n gradient reference value forming portion 22A. 

In this case, the vehicle speed V input to the slip speed calculating portion 20A corresponds to a value calculated 
by the vehicle speed calculating portion 18A. Concretely speaking, it is calculated in the following manner. The slip 
speed Av calculated by the slip speed calculating portion 20A mentioned above and the wheel speed v<d detected by 
the wheel speed sensor 28 are input to the vehicle speed calculating portion 18 A. The sip speed Av corresponds to a 
55 value (V - Vco) obtained by subtracting the wheel speed V© from the vehicle speed V. Accordingly, the vehicle speed 
calculating portion 18A calculates the vehicle speed V by adding the input slip speed Av and the input wheel speed V© 
and outputs it to the slip speed calculating portion 20 A. 

Further, it is known that the road surface \i and the braking force P have an inherent relation in accordance with the 
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vehicle speed as shown in Figs. 27 to 29 in correspondence to each of the road surface conditions. 

Then, each of the gradient reference value forming portions 22A of the reference value forming and storing cir- 
cuits 12D, 12S and 121 stores a relational formula between the braking force P and the road surface \i shown in Figs. 
27 to 29 in accordance with the vehicle speed. Accordingly, when the braking force P c and the vehicle speed V are 
s input, the \i gradient reference value forming portion 22A calculates the road surface u corresponding to the braking 
force P c and the vehicle speed V on the basis of the relational formula mentioned above, the braking force P c and the 
vehicle speed V. 

As mentioned above, the slip speed Av (two points, Av 1 , Av 2 ) corresponding to the braking force P c and the vehicle 
speed V is input to the \i gradient reference value forming portion 22A from the slip speed calculating portion 20A, and 

10 the ii gradient reference value forming portion 22A calculates the road surface \i (two points, ^ , corresponding to 
the braking force P c and the vehicle speed V. Then, the \i gradient reference value forming portion 22A calculates the 
reference value of the gradient Gd of the road surface \x by substituting the slip speed Av and the road surface p for the 
formula (16). The wheel speed Vco and the braking force are respectively sampled at every predetermined sampling 
time from the time of starting the brake to the time when a predetermined time elapses by the wheel speed sensor 28 

is and the braking force detecting portion 34, and input to the u gradient reference value forming portion 22A. Then, the 
n gradient reference value forming portion 22A calculates the reference value of the gradient of the road surface ji 
at every time when the braking force changes by a predetermined amount from the time of starting the brake to the time 
when the braking force is applied by a predetermined amount and outputs the time series data of the reference value 
of the gradient of each of the road surface \i to the reference value time series data storing portion 24A. Then, the 

20 reference value time series data storing portion 24A stores the time series data of the reference value of the gradient 
Gd of the road surface \l 

Further, in the resonance gain time series data storing portion 25, the time series data of the resonance gain is 
stored synchronously with the output timing of the time series data of the reference value from the ji gradient reference 
value forming portion 22A to the reference value time series data storing portion 24A. Accordingly, the resonance gain 
25 GdT and the braking force P c sampled at every predetermined sampling time from the time of staring the brake to the 
time when the predetermined time elapses are input to the resonance gain time series data storing portion 25, and the 
resonance gain time series data storing portion 25 stores the resonance gain QjT at every time when the braking force 
changes by the predetermined amount from the time of starting the brake to the time when the braking force is applied 
by the predetermined amount 

da Next an identification operating routine which the time series data track identifying portion 1 4A performs at a time 
of staring the brake will be described below with reference to Fig. 30. 

In a step 82 of Fig. 30, a variable k expressing a time when the braking force changes by a predetermined amount 
from the time of starting the brake is initialized, and in a step 84, one increment is applied to the variable k. 

In a step 86, the gradient G^T(k) of the road surface \i (the resonance gain) calculated by the resonance gain cal- 
35 culating portion 36 at the time expressed by the variable k (k • x, x: a sampling period) is taken in. 

In a step 88, a variable s representing a number of the road surface conditions to be identified is initialized, and in 
a step 90, one increment is applied to the variable s. 

In a step 92, the reference value G^s(k) at the time k among the time series data of the gradient Qd of the road sur- 
face \i stored in the reference value time series data storing portion 24A corresponding to the road surface condition s 
40 is taken in. 

In a step 94, an identifying value Zq is calculated by the formula (17). 

Z 8 =Z 8+ (G d T(k)-G d s(k)) 2 (17) 

45 

In a step 96, whether or not the variable s is equal to or more than a predetermined total number Sq (3 in this 
embodiment) of the road surface conditions is judged. In the case that the relation of the variable s £ the total number 
Sq is not established, since there is a road surface condition for which the identifying value Zg is not calculated at the 
time k, the step goes back to the step 90 and the processes mentioned above (the steps 90 to 96) are performed again. 

so In the case that the relation of the variable s £ the total number Sq is established, since the identifying value Zg at the 
time k for all the road surface conditions is calculated, in a step 98, whether or not the variable k is equal to or more than 
n expressing a time at which a predetermined time elapses from the time of starting the brake is judged. In the case 
that the relation of the variable k £ n is not established, since the predetermined time does not elapse from the time of 
staring the brake, the step goes back to the step 84 and the processes mentioned above (the steps 84 to 98) are per- 

55 formed again. In the case that the relation of the variable k ^ n is established, since the predetermined time elapses 
from the time of starting the brake, in a step 100, it is estimated that the road surface condition corresponding to the 
minimum identifying value z s among the identifying values z 6 (z-, to Z3) is the road surface condition in which the vehicle 
currently runs. 
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Accordingly, the value Zg can be expressed by the following formula (18). 

n 

Z s =£(G d T(k^G d s(k)) 2 (18) 

5 K=1 

As mentioned above, the value Zg is calculated on the basis of the time series data from the time of starting the 
brake to the time when the predetermined time elapses, because the gradient QjT(k) of the road surface \i calculated 
by the resonance gain calculating portion 36 does not always coincide with the gradient G^ s (k) of the road surface \i in 
w the road surface condition in which the vehicle currently runs (in Fig. 31, it is the road surface condition of Snow), as 
shown in Fig. 31 , and it disperses around the reference value of the gradient of the road surface u in the road sur- 
face condition in which the vehicle currently runs. Accordingly, the road surface condition is accurately estimated by 
using the time series data from the time of starting the brake to the time when the braking force is applied by the pre- 
determined amount. 

75 Then, a peak value m^oc of a coefficient of friction corresponding to the estimated road surface condition is output 
to an ABS control portion or the like (not shown) among the peak values jimax of the coefficient of friction which are 
previously stored in correspondence to each of the road surface conditions, and the data expressing the estimated road 
surface condition is output to the slip speed selecting portion 16 A. 

The ABS control portion provided with the peak value umax the time series data track identifying portion 14A 

20 sets the slip ratio giving the input peak value Umax to be a reference slip ratio and controls the braking force in such a 
manner as to follow the peak \t. 

The slip speed selecting portion 16A provided with the data expressing the road surface condition by the time 
series data track identifying portion 1 4A outputs the slip speed Av corresponding to the estimated road surface condi- 
tion among the slip speeds output from each of the slip speed calculating portions 20A of the reference value forming 

25 and storing circuits 12D, 12S and 121 on the basis of the data expressing the input road surface condition. 

As mentioned above, in accordance with this embodiment, since the road surface condition and the peak p value 
of the road surface on which the vehicle currently runs can be accurately calculated, it is possible to perform the vehicle 
stabilizing control (the change of the control gain, the set of control target value) such as VSC, ABS, TRC and the like, 
to warn the driver of the road surface condition and to estimate each of the vehicle conditions such as a vehicle body 

30 lateral slip, a yaw rate and the like. 

In the embodiments mentioned above, although the relational formula between the braking force P and the slip 
speed Av shown in Figs. 24 to 26 and the relational formula between the braking force P and the road surface \i shown 
in Figs. 27 to 29 are stored, the present invention is not limited to this. Since the braking force P c and the gradient of the 
road surface \i have the relation shown in Fig. 32, it may be structured such as to store the relational formula shown in 

35 Fig. 32 in correspondence to the vehicle speed V, and calculate the gradient of the road surface u corresponding to the 
vehicle speed V and the braking force P c from the vehicle speed V, the braking force P c and the relational formula 
between the braking force P c and the gradient of the road surface u~ 

Further, in the embodiment mentioned above, although the road surface condition is estimated by the minimum 
value of a square error of the gradient (VT(k) of the road surface \i and the reference value O^k), the present inven- 

40 tion is not limited to this. It may be structured such that the road surface condition is estimated by calculating the corre- 
lation between the time series data of the gradient QjT(k) of the road surface \i and that of the reference value GdS(k) 
from the time of starting the brake to the time when the predetermined time elapses. 

Further, in the embodiment mentioned above, the structure is made such as to store the relational formula between 
the braking force P and the slip speed Av shown in Figs. 24 to 26, the relational formula between the braking force P 

45 and the road surface \i shown in Figs. 27 to 29, or the relational formula between the braking force P and the gradient 
Qd of the road surface n shown in Fig. 32. However, the present invention is not limited to these. Since each of the rela- 
tions between the slip speed Av and the road surface \i and between the slip speed Av and the gradient Gd of the road 
surface \i becomes similar to those between the braking force P and the slip speed Av, the road surface \i, or the gradi- 
ent Qd of the road surface \u the structure may be made such as to store the relational formulas between the slip speed 

so Av and the road surface \i and between the slip speed Av and the gradient of the road surface \i and process in the 
same manner. 

Further, in the embodiments mentioned above, the resonance gain is calculated as the gradient of the road sur- 
face it. However, the present invention is not limited to this. Since the resonance gain is a physical amount equivalent 
to the braking torque gradient, the structure may be made such as to include the braking torque gradient calculating por- 
55 tion in place of the resonance gain calculating portion 36, to calculate the braking torque gradient in place of the gradi- 
ent of the road surface u, and to process in the same manner. 

Next a method of calculating the braking torque gradient will be described below. 

A wheel motion of each of the wheels and a vehicle body motion can be described by an equation of motion shown 
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by the following formulas (19) and (20). 

J<s>i = R c F i , (v/R c -© i )-T bi (19) 

5 4 

M v = "S F J , ( v/R c-«() (2°) 

10 In which Fj' is a braking force generated in a wheel of number i, is a braking torque applied to the wheel of 
number i in correspondence to a stepping force, M is a vehicle mass. Rc is an effective radius of the wheel, J is a wheel 
inertia and v is a vehicle speed. In this case, a mark • denotes a differential with respect to a time. In the formulas (19) 
and (20), Ff is shown as a function of the slip speed (v/Rc - ©j). 

Here, the vehicle speed is equivalentty expressed by an angular velocity of the vehicle body (the formula (21)), 

75 and the braking torque R^' is described as a linear function of the slip speed (gradient kf, y component T) (the formula 
(22)). 



"C V 

20 



v=R c g> v (21) 
R c F|^-»i>-^ W®!^, (22) 



Further, the equations (23) and (24) are obtained by substituting the formulas (21) and (22) for the formulas (19) 
25 and (20) and to express the wheel speed cDj and the vehicle speed av as the time series data ©fk] and ovfkj (krsa sam- 
ple time in the case of setting the sample time x to be a unit, k = 1 , 2, ... ) dispersed at every sample times x. 

cj>j[k]-G>j[k-1] 

J ' / = k^^k-IJ-cDilk-llHTrTy (23) 

30 

R c 2 M G>v[kl 7 vlk " 11 - -S k i ' •v^1]+Z«k l «,lk-1])-£T J (24) 
H H H 

35 

Here, when the formulas (23) and (24) are made simultaneous and the equivalent angular velocity ©y of the vehicle 
body is deleted, the following formula (25) is obtained. 

4 4 

R c M J=1 R c M i=1 

2 4 2 4 2 4 ^ 

'jr c 2 m~ J J jr c *m£ j jr c 2 m^ j 

45 

In this case, when it is supposed that the maximum braking torque of RcMg/4 (g is a gravitational acceleration) is 
generated under the condition of the slip speed 3 rad/s, the following formula (26) is obtained. 

so R c Mg/4 

max (k,)- (26) 



In this case, as a concrete constant, when x = 0.01 (sec), Rc = 0.3 (m) and M = 1000 (kg) are considered, the fol- 
55 lowing formula (27) is obtained. 
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4 



rnax(— ^^* y ) * 0.054 « 1 (27) 



5 



10 



15 



Accordingly, the formula (25) can be approximated as the following formula (28). 

k,j{<D,[k-1]-a>,Ik-2]}+f i = -o,[k]-i2cD,[k-1l-G),[k-2] (28) 

In which f j is expressed by the following formula (29). 



In the case that a brake hydraulic pressure (a master cylinder pressure or a wheel cylinder pressure) can be meas- 
20 ured with respect to T b in the formula (29), the following formula (30) can be established. 

T b = k b -P c (k) (30) 

25 In which P c is a brake hydraulic pressure and k^ is a coefficient for converting the brake hydraulic pressure into the 
braking torque. In this case, the formulas (28) and (29) are converted into the formulas (31) and (32). 

kij{G>i[k-1]-Oi[k-2]}+f , = -a> i [k}+2a) j [k-1]-<D i [k-2]-P ci [k-1] (31) 

30 

2 4 2 4 

JR C 2 M~ J JR C 2 M~ J 

35 

In this case, P d is obtained by the following formula (33). 

2 4 

Pc,[k] = — L_2kj.T b (h) (33) 
40 JR C 



Here, it is supposed that the brake hydraulic pressure P c (k) at the present sampling time k and the brake hydraulic 
pressure P c (k-1) at a time k-1 of one sampling prior thereto are substantially equal to each other (p c (k) = P c (k-1) ). 
45 Further, the preceding estimated value is used as the ji gradient kj in the formula (33). 

The formula (31) can be described in a form of a linear manner with respect to unknown quantities kj and fj by 
arranging in this manner, and the braking torque gradient kj with respect to the slip speed can be estimated by applying 
the on-line parameter identifying method to the formula (31). 

Accordingly, the time series data of the braking torque gradient can be estimated from the time series data <Dj[k] of 
so the detected wheel speed by repeating the following steps 1 and 2 (method of least squares). 

Stepl: 

*,M T -e f = y| [k] (34) 
In which the following formulas (35), (36) and (37) are established. 
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35 



40 
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15 yjM^-^ilk^cDjlk-ll-cajlk^l-PcjCk) (37) 

In this case, a first element of the matrix <fc[k] in the formula (35) is a physical amount with respect to the change of 
the vehicle speed during a sample time, and the formula (37) is a physical amount with respect to the change in the 
20 change of the wheel speed during the sample time (a physical amount with respect to a value obtained by differentiating 
the wheel speed two stages). 

Step 2: 

25 QJk] = eilk-IJ+L^KyilkHiM 7 ^^-!]) (38) 



^+*iM T P i tk-1]* i M 



L|M - _'t_ . .. - ( 39 > 



4 *+W} T W-Mlk} J 



(40) 



An estimated value of 0j is calculated by the above recurrence formulas, so that a first element of the matrix of the 
estimated value of 6j is extracted as an gradient of the estimated braking torque. In this case, X is a forgetting factor (for 
example, X = 0.98) showing a degree for removing the past data, and " T " shows a transposition of the matrix. 

In this case, a left side of the formula (38) is a physical amount for expressing a history of the physical amount with 
45 respect to the change of the wheel speed and a history of the physical amount with respect to the change in the change 
of the wheel speed. 

As shown in Figs. 33 to 35, in each of the road surface conditions comprising Dry, Snow and Ice, the time series 
data QdTj of the gradient Gd of the road surface \i calculated by multiplying the braking torque gradient calculated by 
using the formula (29) by a conversion factor accurately coincides with the time series data Gd 8 of the reference value 

so of the gradient Gd of the road surface and is further close to the track of the actual running road surface, so that an 
estimating accuracy of the road surface condition can be improved. Here, Figs. 33 to 35 show the time series data QjT 2 
of the gradient Gd of the road surface \i calculated by multiplying the braking torque gradient calculated by the formula 
(29) by the conversion factor. The time series data GdT-i more accurately coincides with the time series data Gd& than 
the time series data GdT 2 . Accordingly, the time series data GyTj is improved compared to the time series data GdT 2 . 

55 Next a fifth embodiment in accordance with the present invention will be described in detail with reference to the 
drawings. 

As shown in Figs. 39A and 39B,a characteristic between the gradient of the road surface \i and the slip speed has 
a great difference at each of the road surfaces at a time of starting the brake at which the slip speed is small. 
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Accordingly, in accordance with the fourth embodiment, there is suggested a road surface friction condition esti- 
mating apparatus for comparing the gradient of the road surface \i detected in an on-line manner by using the wheel 
speed, the braking force and the like with the reference value of the \i gradient of each of the road surfaces comprising 
Dry, Snow and Ice previously prepared on the plane constituted by the braking force P c and the u gradient so as to esti- 
5 mate the road surface friction condition from an initial time of starting the brake. 

However, there is a relation as shown in Fig. 40 between the slip speed Av and the coefficient \i of friction. In other 
words, there is a time difference between the change of the slip speed and the change of the coefficient \i of friction in 
accordance with an increase of the braking force, so that a micro hysteresis loop is generated. Accordingly, the \i gra- 
dient detected in an on-line manner appears as a great noise component with respect to the braking force in a horizon- 
10 taJ axis as shown in Fig. 41 . 

As a method of removing such a noise at every fixed sampling time At, there is widely known a low pass filter of the 
type shown by the formula (41). 

15 xfi+1) = (1-^)X(i)+Yu(i) (41) 



In this case, X is an output state amount after filter-treatment u is an input state amount on which a filter-treatment 
is performed, T is a time constant for giving a desired filter effect and i is a sampling period renewed at every At. 
20 Fig. 42 expresses a result of shows that fitter-treated in accordance with the formula (41 ), which shows that a noise 
is not sufficiently removed. This is because the filter-treatment is performed with respect to the time change, but without 
relation to the braking force corresponding to the horizontal axis. Accordingly, in comparing with the \i gradient refer- 
ence value, the degree of comparison is lowered, so that there is a risk of deteriorating an accuracy of estimating the 
road surface condition. 

25 This embodiment is made by taking the fact mentioned above into consideration, and is structured such as to accu- 
rately estimate the road surface condition. 

This embodiment will be described in more detail below. In the road surface condition estimating apparatus in 
accordance with this embodiment (to which a variation reduction operating apparatus is applied) has the same compo- 
nents as those of the road surface condition estimating apparatus in accordance with the fourth embodiment mentioned 

30 above, so that the same reference numerals are attached to the same elements and the description thereof will be omit- 
ted. 

The road surface condition estimating apparatus in accordance with this embodiment is different from the fourth 
embodiment mentioned above in view that as shown in Fig. 36, the resonance gain storing circuit 12j is provided with 
a resonance gain fitter portion 26A for performing a fitter-treatment at every time when the braking force changes by a 

35 predetermined amount with respect to the calculated resonance gain, and a resonance gain time series data storing 
portion 25A for storing a time series data of the resonance gain filter-treated by the resonance gain filter portion 26A. 
The braking force detecting portion 34 and the resonance gain calculating portion 36 are connected to the resonance 
gain filter portion 26A. The resonance gain filter portion 26A is connected to the resonance gain time series data storing 
portion 25A. The resonance gain time series data storing portion 25A is connected to the time series data track kJerrti- 

40 tying portion 14A. 

Next, an operation of this embodiment will be descrfoed below. 

Since a method of forming a reference value of the gradient of the road surface u by the reference value forming 
and storing circuits 1 2D, 1 2S and 1 2 1 in accordance with this embodiment is the same as the forming method in accord- 
ance with the fourth embodiment mentioned above, the description thereof will be omitted. 

45 In the resonance gain time series data storing portion 25A, the time series data of the resonance gain is stored in 
synchronously with the output timing of the time series data of the reference value to the reference value time series 
data storing portion 24A from the ji gradient reference value forming portion 22 A. Accordingly, the resonance gain G^T 
fitter-treated in the resonance gain filter portion 26A is input to the resonance gain time series data storing portion 25A 
at every time when the braking force changes by a predetermined amount from the time of starting the brake to the time 

so when the braking force is applied by the predetermined amount mentioned above, and the input resonance gain GyT is 
stored therein. 

A problem of the fourth embodiment mentioned above exists in a point that a fitter-treatment is performed at every 
sampling period. Then, a dispersion of the track in Fig. 42 becomes small by performing the fitter-treatment only when 
the braking force changes by a predetermined amount. Concretely speaking, the filter-treatment expressed by the fof- 
55 lowing formula (42) which is similar to the formula (41 ) is performed. 

X( k+1) = (1 -^5) - X (k)+^u(k) (42) 
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In which Tp is a constant corresponding to the time constant T in the formula (42), AP C is a braking force of a pre- 
determined amount and k is a period renewed at every time when the predetermined braking force changes by a pre- 
determined amount. The equation may be expanded so as to be expressed by the following formula (43). 
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In which X is a variable vector of an internal condition of the filter (size is n x 1), and a1 1 to arm, b1 to bn, c1 1 to 
c1n and d1 to dn are respective constant of the filters. In Fig. 43, the filter effect of the formula (42) is compared with 
that of the formula (41). As shown by Fig. 43, it is understood that the dispersion is restricted to a smaller amount in 
comparison with the filter-treatment result in accordance with the formula (41), and an average value as a result of the 

20 formula (41) can be obtained at each of the points of the braking force. In this case, at a time of starting the braking 
force (adjacent to the braking force 0), it moves along the maximum value of the gradient of the road surface |i. This 
shows an excessive response condition with respect to the magnitude of Tp. 

In place of the fitter-treatment as shown by the formulas (42) and (43), an average value during the period in which 
the braking force changes by a predetermined amount may be calculated in a manner shown by the following formula 

25 (44). 

m 

m - E(u©ym (44) 

30 

In which m is a number of resonance gain points sampled from k to k +1 , x(k) is an average value of a resonance 
gain at a time of k, and j is a period initialized to 1 at a time of k and renewed at every sampling periods. 

Next, the identification operating routine performed by the time series data track identifying portion 1 4A at a time of 

35 starting the brake will be described below. The identification operating routine in accordance with this embocfi merit is 
the same as the identification operating routine in accordance with the fourth embodiment (refer to Fig. 30). However, 
in a step 86, the gradient GdT(k) of the road surface \i (the resonance gain) stored in the resonance gain time series 
data storing portion 25A at a time expressed by the variable k (k-x, x: optional multiple of a sampling period) is taken 
in. The gradient QdT(k) of the road surface u is obtained by f Bter-treating on the basis of the formula (42) or (43) by the 

40 resonance gain fitter portion 26A. 

As mentioned above, in accordance with this embodiment, since it is possible to accurately calculate the condition 
of the road surface on which the vehicle currently runs and the peak \i value, it is possible to perform a stabilizing control 
of the vehicle (change of the control gain and set of the control target value) such as VSC, ABS, TRC and the like, to 
warn the driver of the road surface condition and to estimate each of the vehicle conditions such as the vehicle lateral 

45 slip angle, the yaw rate and the like. 

In the embodiment mentioned above, the resonance gain is filter-treated by using the formulas (42) or (43). Alter- 
natively, the structure may be made such as to calculate an average of the resonance gains C^T output at every sam- 
pling periods from the variable k to k + 1 . 

Further, in accordance with this embodiment since the braking force P shown in Hgs. 24 to 26 and the gradient of 

so the road surface \i have the relation shown in Fig. 32 in the same manner of the fourth embodiment mentioned above, 
the structure may be made such as to store the relational formula shown in Fig. 32 in correspondence to the vehicle 
speed V and calculate the gradient of the road surface n correspond rig to the vehicle speed V and the braking force P c 
in accordance with the vehicle speed V, the braking force P c and the relational formula between the braking force P c 
and the gradient of the road surface u- 

55 Still further, in accordance with the embotfment mentioned above, the road surface condition is estimated in 
accordance with the minimum value of the square error between the gradient G^T(k) of the road surface ji and the ref- 
erence value G^k). However, the present invention is not limited to this. The structure may be made such as to calcu- 
late the correlation between the time series data of the gradient G^T(k) of the road surface ji and that of the reference 
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value G d3 (k) from the time of starting the brake to the time when the predetermined time elapses so as to estimate the 
road surface condition. 

Furthermore, in accordance with the embodiment mentioned above, there are stored the relational formula 
between the braking force P and the slip speed Av shown in Figs. 24 to 26 , the relational formula between the braking 

5 force P and the road surface n shown in Figs. 27 to 29 , or the relational formula between the braking force P and the 
gradient of the road surface \i shown in Fig. 32. However, the present invention is not limited to these. Since each of 
the relations between the slip speed Av and the road surface u and between the slip speed Av and the gradient Qd of 
the road surface n have similar to those between the braking force P and the slip speed Av, the road surface ji, or the 
gradient of the road surface \i, the structure may be made such as to store the relational formulas between the slip 

10 speed Av and the slip speed Av, and the road surface \i and between the slip speed Av and the gradient of the road 
surface n, and to process in the same manner. 

Moreover, in the embodiment mentioned above, the resonance gain is calculated as the gradient Gd of the road 
surface \i. However, the present invention is not limited to this. Since the resonance gain is a physical amount equivalent 
to the braking torque gradient, the braking torque gradient calculating portion is provided in place of the resonance gain 

is calculating portion 36, the braking torque gradient is calculated in place of the gradient of the road surface u, and to 
process in the same manner. 

As shown in Fig. 43, in the track of the resonance gain filter-treated in accordance with the formula (42), it is under- 
stood that the dispersion is restricted to be smaller at each of the braking force points than the track filter-treated in 
accordance with the formula (41). 

20 Then, as shown in Figs. 37 and 38, in each of the road surface conditions comprising Dry and Snow, the time series 
data Q d T 1 of the gradient of the road surface u calculated by multiplying the braking torque gradient calculated in 
the manner mentioned above by the conversion factor and filter-treated in accordance with the formula (42) accurately 
coincides with the time series data of the reference value of the gradient Gd of the road surface n, and is further 
close to the track of the actual road surface on which the vehicle runs, so that an estimating accuracy of the road surface 

25 on which the vehicle runs, condition can be improved. 

In accordance with the embodiment mentioned above, the braking force to the wheel is used as the first physical 
amount for expressing the braking condition or the driving condition of the wheel in the vehicle. However, the present 
invention is not limited to this, and the sip speed can be employed. 

Further, in the embodiment mentioned above, the variation reduction operating apparatus is applied to the road 

30 surface condition estimating apparatus. However, the present invention is not limited to this, and the variation reduction 
operating apparatus can be also applied to the ABS apparatus. 

Claims 

35 1. A road surface condition estimating apparatus provided with a correlation among a plurality of wheel motion 
amounts expressing the characteristic of a coefficient of friction \i between a wheel and a road surface with respect 
to a slip speed for each road surface condition, comprising: 

memory means for storing a gradient of the coefficient of friction u with respect to the slip speed or a physical 
40 amount relating to said gradient as one of the wheel motion amounts expressing said correlation; and 

road surface condition estimating means for estimating and oulputting a road surface condition on the basis of 
comparison between a converted value of the wheel motion amount for each road surface condition obtained 
by respectively converting at least one of detected values of the plurality of wheel motion amounts in accord- 
ance with said stored correlation and a detected value of the same amount as said wheel motion amount 

45 

2. A road surface condition estimating apparatus comprising: 

detecting means for detecting a plurality of wheel motion amounts expressing a characteristic of a coefficient 
of friction u between a wheel and a road surface with respect to a slip speed; 
so converting means for converting a detected value of at least one of the plurality of wheel motion amounts into 

a converted value of another of the plurality of wheel motion amounts for each of a plurality of road surface con- 
ditions on the basis of a correlation between the plurality of wheel motion amounts for each of a plurality of road 
surface conditions provided in said converting means; and 

road surface condition estimating means for estimating a road surface condition by comparing said converted 
55 value with a detected value of said another of the plurality of wheel motion amounts. 

3. A road surface condition estimating apparatus according to claim 2, wherein one of said plurality of wheel motion 
amounts is a physical amount expressing slipperiness of the wheel. 
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4. A road surface condition estimating apparatus according to claim 3, wherein said physical amount is a gradient of 
a coefficient of friction \i with respect to the slip speed or a physical amount relating to said gradient 

5. A road surface condition estimating apparatus according to claim 4, wherein said physical amount relating to said 
5 gradient is a resonance gain. 

6. A road surface condition estimating apparatus according to claim 2, further comprising slip speed calculating 
means and which speed calculating means wherein 

w said detecting means detects a wheel speed as one of said plurality of wheel motion amounts, 

said correlation includes a correlation defined for each vehicle speed and a correlation between two wheel 
motion amounts other than the vehicle speed, 

said converting means converts a detected value of any one of said two wheel motion amounts into a con- 
verted value of the other wheel motion amount for each of a plurality of slip speeds on the basis of said corre- 
15 lation between said two wheel motion amounts, 

said slip speed calculating means calculates a slip speed on the basis of said converted value and a detected 
value of said the other wheel motion amount, 

said vehicle speed calculating means calculates the vehicle speed on the basis of the slip speed calculated by 
said slip speed calculating means and the wheel speed detected by said detecting means, and 
20 said road surface condition estimating means estimates the road surface condition on the basis of the vehicle 

speed calculated by said vehicle speed calculating means and said correlation defined for each vehicle speed. 

7. A road surface condition estimating apparatus according to claim 2, further comprising peak u value searching 
means for searching a peak \i value corresponding to said estimated road surface condition on the basis of peak \i 

25 values provided for each of said plurality of road surface conditions and said road surface condition estimated by 
said road surface condition estimating means. 

8. A road surface condition estimating apparatus according to claim 2, wherein a maximum braking force is provided 
for each of a plurality of road surface conditions and a plurality of vehicle speeds, and 

30 

said detecting means detects at least vehicle speed and a braking force as said plurality of wheel motion 
amounts, and wherein said road surface condition estimating apparatus further comprises: 
maximum braking force selecting means for selecting the maximum braking force corresponding to said esti- 
mated road surface condition and said detected vehicle speed on the basis of the maximum braking force pro- 
as vided for each of said plurality of road surface conditions and said plurality of vehicle speeds, the road surface 
condition estimated by said road surface condition estimating means and the vehicle speed detected by said 
detecting means, and 

subtracting means for calculating the difference between the maximum braking force selected by said maxi- 
mum braking force selecting means and the braking force detected by said detecting means. 

40 

9. A road surface condition estimating apparatus according to claim 6, wherein a maximum braking force is provided 
for each of a plurality of road surface conditions and a plurality of vehicle speeds, and 

said detecting means detects the braking force as one of said plurality of wheel motion amounts, and wherein 

45 said road surface condition estimating apparatus further comprises: 

maximum braking force selecting means for selecting the maximum braking force corresponding to said esti- 
mated road surface condition and said calculated vehicle speed on the basis of the maximum braking force 
provided for each of said plurality of road surface conditions and said plurality of vehicle speeds, the road sur- 
face condition estimated by said road surface condition estimating means and the vehicle speed calculated by 

so said vehicle speed calculating means, and 

subtracting means for calculating the difference between the maximum braking force selected by said maxi- 
mum braking force selecting means and the braking force detected by said detecting means. 

10. A road surface condition estimating apparatus comprising: 

55 

first detecting means for detecting a first physical amount expressing a braking condition of a wheel of a vehi- 
cle; 

second detecting means for detecting a vehicle speed of said vehicle; 
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third detecting means for detecting a second physical amount which represents slipperiness of said wheel and 
which is determined in accordance with a road surface condition on which said wheel is running, said first phys- 
ical amount, and said vehicle speed; 

calculating means for calculating said second physical amount for each of a plurality of predetermined road 
5 surface conditions on the basis of said first physical amount and said vehicle speed; 

comparing means for comparing a detected track of said second physical amount detected by said third detect- 
ing means with respect to said first physical amount with each of calculated tracks of said second physical 
amounts calculated for each of the plurality of road surface conditions by said calculating means with respect 
to said first physical amount; and 
10 estimating means for estimating a road surface condition corresponding to a track closest to said detected 

track among said calculated tracks as an actual road surface condition in which said wheel is running on the 
basis of the comparison results of said comparing means. 

1 1 . A road surface condition estimating apparatus according to claim 1 0, wherein said first physical amount is a braking 
is force on the wheel or a slip speed. 

12. A road surface condition estimating apparatus according to claim 10, wherein said second physical amount is any 
one of a gradient of a coefficient of friction u between the road surface and the wheel with respect to the slip speed, 
a gradient of braking torque on the wheel with respect to the slip speed, and a gradient of braking force on the 

20 wheel with respect to the slip speed. 

13. A road surface condition estimating apparatus according to claim 10, wherein 

said first detecting means, said second detecting means and said third detecting means respectively detect the 
25 first physical amount, the vehicle speed, and the second physical amount at each predetermined time from the 

time braking is started on a wheel, 

said calculating means calculates the second physical amount for each of a plurality of road surface conditions 
at each time the first physical amount changes by a predetermined amount from the time braking is started on 
the wheel, and 

30 said comparing means compares the second physical amount detected at said predetermined time, with each 

of the second physical amounts calculated for each of the plurality of road surface conditions and at each time 
the first physical amount changes by a predetermined amount at each time the first physical amount changes 
by a predetermined amount from the time braking is started on a wheel. 

35 14. A variation reduction processing apparatus comprising: 

first detecting means for detecting a first physical amount expressing a braking condition or a driving condition 
of a wheel of a vehicle; 

second detecting means for detecting a second physical amount which represents slipperiness of said wheel 
40 and varies in accordance with a change in said first physical amount; and 

reduction processing means for reducing variations in the second physical amount in accordance with said 
change in the first physical amount. 

15. A variation reduction processing apparatus according to claim 14, wherein said reduction processing means per- 
45 forms a filter-treatment on variations in said second physical amount at each time said first physical amount 

changes by a predetermined amount, or calculates an average of said second physical amount during the time said 
first physical amount changes by a predetermined amount 

16. A road surface condition estimating apparatus comprising: 

so 

first detecting means for detecting a first physical amount expressing a brake condition or a driving condition 
of a wheel of a vehicle; 

second detecting means for detecting a vehicle speed of said vehicle; 

third detecting means for detecting a second physical amount which represents slipperiness and which is 
55 determined in accordance with a road surface condition on which said wheel is running, said first physical 

amount, and said vehicle; 

reduction processing means for reducing variations in the second physical amount detected by said third 
detecting means in accordance with a change in said first physical amount; 
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calculating means for calculating the second 
physical amount for each of a plurality of prede- 
termined road surface conditions on the basis 
of the first physical amount and the vehicle 
speed; 

comparing means for comparing a detected 
track of the second physical amount in which 
variation is reduced by said reduction process- 
ing means with respect to the first physical 
amount with each of calculated tracks of the 
second physical amount calculated for each of 
the plurality of road surface conditions by said 
calculating means with respect to said first 
physical amount; and 

estimating means for estimating a road surface 
condition corresponding to a track closest to 
said detected track among said calculated 
tracks as an actual road surface condition on 
which said wheel is running on the basts of the 
comparison result of said comparing means. 

17. A road surface condition estimating apparatus 
according to claim 16, wherein said reduction 
processing means performs a filter-treatment on 
variations in said second physical amount at each 
time said first physical amount changes by a prede- 
termined amount or calculates an average of said 
second physical amount during the time said first 
physical amount changes by a predetermined 
amount 

18. A road surface condition estimating apparatus 
according to claim 16, wherein said first physical 
amount is a braking force on the wheel or a slip 
speed. 

19. A road surface condition estimating apparatus 
according to claim 16, wherein said second physi- 
cal amount is any one of a gradient of a coefficient 
of friction between the road surface and the wheel 
with respect to the slip speed, a gradient of braking 
torque on the wheel with respect to the slip speed, 
and a gradient of braking force on the wheel with 
respect to the slip speed. 

20. A road surface condition estimating apparatus 
according to claim 16, wherein 



amount from the time at which braking is 
started on the wheel, and 
said comparing means compares the second 
physical amount detected at said predeter- 
mined time with each of the second physical 
amounts calculated at each time the first phys- 
ical amount changes by a predetermined 
amount, at each time the first physical amount 
changes by a predetermined amount from the 
time at which braking is started on the wheel. 



said first detecting means, said second detect- 
ing means and said third detecting means 
respectively detect the first physical amount 
the vehicle speed, and the second physical 
amount at each predetermined time from the 
time at which braking is started on the wheel, 
said calculating means calculates the second 
physical amount for each of a plurality of road 
surface conditions at each time the first physi- 
cal amount changes by a predetermined 
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